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SHOCK TUNNEL STUDIES OF SCRAMJET PHENOMENA 1993 
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I mA with the format of previous reports, this consists of a series of reports on 
specillc projects. There is a brief introduction commenting on each report, and the project 
reports follow in the order of the headings in the introduction. 


-Program A" corresponds to work funded jointly by NAGW-674 and Australian sources^ 
and "Program 13" to that funded by Australian sources alone. 


PROGRAM A 


(0 


COMPARATIVE S T UDIES OF T4 jVN PJIYPlJ jLSE 


Chemical Recon, hi nation in an E xpansion Tube ( R. Bahos a„d R.G. Morgan) 

... . ;„„i cit.Hu hac been carried out to determine the probable "frozen" 

rs? a ror diap,,ra6m rup,u vr 

■ j u -j nB this model it was found that the test gas composition could be 
obtained from an entropy correla'.ion, similar to that employed by Harris for stead, nozzle 

flows. 


(ii) DIRECT MEASUREMENT OF T11RJJ ST/DRAC 


(a) Drag Meas urejTLent^m-muntedX^n es and a Scramj gL^cJiLcje^ lyperyelodty 
plow (L. Porter, D. Mee and A. Pauli) 

This records our first measurement of drag on a scramjet model (as well as drag 
measurements on blunt cones). Drag coefficients of the model were found to be qui e 
high, and were subsequently reduced by reducing the afterbody cone angle, and the ang 
of the afterbody of die splitters between the combustion chambers. 

(b) Tbilisi Measurements of a Co mEki^Ajiiammciri 

(A. Pauli, R.J. Stalker and D.J. Mee) 

This records the early measurements of thrust when fuel was injected on the scramjet 
model above. The model was modified by reducing the thrust surface angles as noted. 
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(C) 


cimrW Tunnel (A Pauli, R.J. Stalker and D.J. 
^.rom jpt Thrust Mca surgmenUn.^ '■ 

Mce) 


|llvv ) 

The effect of stagnation ’^gnX e'lalpyfs 

«»< the positive ^ thrust is recorded. Subsequent investigate 

sTrowetHhaUhis^vas primarily due to the increase in precombnstion temperature. 

(d) Thrust Measurerng nUO-g^^^ < S ' DJ ' M " ^ 

Simmons) 

Consider* difHcu.ty has T, 5 

X i» l. *m P friction iosses do no, make a iarge difference to 

ilip thrust of the nozzle. 


(Hi) rOMRUSTlQIl SCALlNG 


Scaling in a sctanri e t j^nbnsior (M.V. Pulsonetti and R.J. Stalker) 

* j' i* „ #uj»* cralinff in scrjiniict combustors may 
rUne- tha, depends on the product 

^ararieristic dimension). 


follow that of 
(pressure x a 


PROGRAMS 


(j) mixing AND i ^MBUSIIQN 

(a) jJypersPBicJgnittfiiLSDdJillllSLEffi^ySfi^lOlk^^^®^ (A ra " n) 

U H,.rt with a nozzle attached have been made. 

Measurements on a long com usi obtained with hydrogen fuel, but this 

Specific impulse values in excess of 1500 sec were omam 

decayed rapidly as stagnation enthalpy was increase . 

(b) Su personic C ombustio n . a Constan LArga_Duct (M. Wendt) 

This began as a PhD study °[ th ^ a substantia^ 

ZuntTf w“k on^otheTeffects was done. Ue work has been presented as a , bests, and 
the conclusions only are reported here. 
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(c) Srramiet Cnmbuslion Surveys (K. Skinner) 


A PhD thesis on the mass spectrometer is currently being examined, and a chapter from 
that thesis relating to mass spectrometric 

« _ - « *. 

profiles, indicating that species and momentum diffusion are similar. 


(ii) SKIN FRICTION AND HEA T TRANSFERMEASUREMENT S 
(a) Tranciiinn nf ComnressiblcJIi g h finfhalpy Boun da t xL^^^^ 

(Y. He and R.G. Morgan) 

, f „ , _ t fran „ fpr to a nat plate in T4 showed transition taking place at 

ReynoldTlIumber ’between 2,8 x 10‘ and 0.6 x 10* for stagnation enthalpies arranging 
from 3 to 7 MJ/kg. 

Skin friction measurements in a turbulent boundary layer are being attempted. 


(iii) EXPANSION tub es 


Sn pprnrhital ExpansionJ ^be^no^ (A- Neely and R.G 

Morgan) 


This paper reports experiments with a small expansion tube in which a usable test flow 
was produced at a test section velocity of 1 3 km/sec. 


(iv) FORCE BALANCE 

Balances for the Measureme nLaLMuhijileC^ 

n^aihm Tb.J. Mee, W.J.T. Daniel, S.L. Tuttle and J.M. Sunmons) 

The principle of the force balance, which has been used to measure thrust and drag in 
shock tunnels, is extended to measure lift and pitelnng moment as well. 
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Chemical Recombination in an 
Expansion Tube 

Robert J. Bakos* and Richard G. Morgant 
University of Queensland. 

Brisbane. Queensland 4072, Australia 

Introduction 

Si^,oi. i. ui 

•ram (Fig la), this requires that the secondary diaph f a8 
roDWre instantaneously on impact by the primary shock, and 
that its mass contributes no inertia to the test gasfor subse- 

bination*in the subsequent unsteady expansion. In this sense, 
the problem is similar to the nozzle freezing problem occurring 

■•asttf srsisssst. ^ - -s 

bo h the strength of the reflected shock generated and the rate 
tfw^nsion experienced by the test gas as it . «cc und " 
■he influence of the unsteady expansion. This Nme •dopis | 
inertia rupture model 3 assuming that the di 
* tfrapm shears cleanly around its periphery, remains intact 
S «,l « 5provid« no resistance other than its 

and nearly planer p ^ preJsufe fie , d The diaphragm 

iner ,a relieves the pressure behind the reflected shock so 
^h« rweaket ”.h «ime (Fi,. lb), lowering the test gas 
enuojy and initial dissociation levels. Also, and perhaps more 
importantly, the diaphragm inertia limits the initial expansion 
taie allowing effective recombination to be achieved. 

F„*rS.'io« „ ,ppro,im,ie 

Laarangian pressure lime history is presented for tn '«***“• 
Aone-dimensional numerical solution with equilibrium chem- 
i“y for the diaphragm-inertia model was also done and shows 
goJd agreement with both the approximate analyncalsolutmn 

£da pressure measurement for air test gas from the NASA 
‘ ~ TT < .oov rrvUion received Dec. 29, 1993; accepted for 

^Senior L^urer. Department of Mechanical Engineering. Member 
AIAA. 


Langley HYPULSE expansion tube. 4 The extent of recombi- 
nation is then calculated using a finite rate chemistry P* c k*g< 
for the analytically or numerically determined pressure-time 

^Calculations including the diaphragm inertia show signifi- 
cantly lower dissociation levels than those reported recently by 
Wilson’ where the diaphragm was taken to hold fast for a 
specified holding lime and then allowed to open instanta- 
neously without contributing inertia to the subsequen accelen 
ation (Fig- 1c). In Ref- 5. this model was used in a full 
Mimrrmonal unsteady-flow calculation with finite rate 

fois rlpVure model b included here for comparison wtth the 

* n FinalW^two-parameter correlation summarizes the re™'* 5 
onhe nonequilibrium calculation, for air in term, of the 




n« I »-/ diagrams lor secondary diaphragm raplarr models: a) Ideal 
diaphragm, b) Inertial diaphragm aad e) boldlag rime diaphragm. 
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josuefkcted of 

ime *cak for the diaphragm. Th» ^ of diaphragm 

lonequilibrium air Comptons <« *"£ * lube . 

and suborbital test conditions m an expa 

Analysis 

Under the assumption of a it** inertia and 
the diaphragm m<^on wiUp the diaphragm 

ihe pressure applied to * uher lhc twl and acceleration 

acts as a very light piston «g»J ra.ios'.re of 

gases. Because typical pr 1 in P (he * ompr „sing acceleration 
the order 100-1000:1. pres (ion 0 f the diaphragm 

•“ b r*Wf p„ •>* «• 

iSSnsST-n. ,1 “ p '"“ ,n 

du ^ Pr (1) 

df Pr w r 

where P* is the dia P’’ r **" ^MI k diaphragm and the 
Following Meyer, 8 as ^ # uniform value of the 

reflected shock is ¥ sum !* / , _ u + u> where a is the sound 
Ricmann invariant J+ - . suc h t ^ at the /- charac- 

speed and u the * 4 *. mot ‘? m analing from the diaphragm 
terisucs will be straight ^“‘"Xough , simple wave, 
surface. That is, the exp either the reflected shock be 

Uniformity of J. requlr “ * aC ross it can be ignored, or 
weak such that entropy change* ^ enttopy region 

that it does not decay so th* **J . t , he initial renected 
exists behind it. For conditions o t within the range 

shock Mach number relative to he t«t 8^ ^ iignincant ,y 

M n - 2-3. which, althoughy^a^ |he (ully reflected 

alter the value of J. • . h , he diaphragm motion is 

shock and assumed not .' 0 ' ' ^ inr.nite length reservoir at the 
equivalent to that cause y «ound speed. This leads to 

the following relation for the preu 
the <Jiaphragm, 

sK^r 

3i‘2K5SSi A <■> - ■* « 

the diaphragm velocity 

a.,. (,.f) 01 

,h, Z-™ «'»> •> ■ p *“ u 

diaphragm face, 

+ *> 


*», o,cu„ to * .« 

2BS®SK£.*XU— -w 

£ = 0^*" (5) 

where po '* ‘he values of to- 

Different particle h'Stones P° s d p/df = - 2 ifh/iy + »'o. 

The initial expansion ra e • j ace(|t l0 lhe diaphragm at 

so that the partrek initially ' . ^ an infinite expan- 

rupture. corresponding to 1,-0, experiences 

sion rate. 

Inertial Diaphragm Model Results 
The nominal Mach 17 

PULSE expansion tube jj t J r t0 3.45 kPa. and the 

lion the shock tube .* , 6?0 approaching the diaphragm, 
primary shock speed filled 10 7 2 Pa, also with air, and the 
The acceleration tube is fit * Equilibrium conditions 

behind the renected sh ^ y^e one-dimensional 

MPa, «,= 1380 ^. ^ 7 Th« h . s ^ ^ usi 
finite volume code of Jac °j> ^ so , ution was started 

approximate ^urlibnumchem^ ^ ^ diaphragm and 

at the time the pr ry assumed to act as a 

continued for 0.5 SnvSence was confirmed by 


continued for 0 5 ms. Th^ajragm was ~ ^ by 

and noting no signincant change in the 

solution. j . w- w al| pressure trace’ at a 

Figure 2a shows a measured . p diaphragm with 

location 76 mm ^‘TZuuSonT^eTumerical solution 
analytical and numerical s imul»tion is. 1 ne u , he 

.«•>■•••>» •“*”*•**: 35- rS* .«*.« » « 

initial pressure field. Tn 1 possibly as a result of 

SS’d'»w« w °nJno3$!jo 

ssiSf ssss;-** ,n F,t a 


2-0 *i)' 


,1 r - 0 is dp^/dr = -iPl/Od Td°a, 
The initial expansion r.«.‘» " "“^gm mass, 
which is seen to location and pressure-time 

For subsequent partwl«.^hw toad 0 ^ ^ di- 

history is evaluated an^ytic y ht charac i„,«,a.’ 

ss-r ifS-sss 
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shows that there is good agreement between the numerical and 
analytical estimates of pressure decay for the particle adjacent 
to the diaphragm. This also indicates that the effects of chem- 
istry on the pressure distribution are small for the conditions 
considered and that the ideal gas approximation is adequate 
for this study. 


Recombination Chemistry 

Recombination in the expansion is calculated using a gen- 
eral chemical kinetics computer program 10 for a specified par- 
ticle pressure-time history, with the gas Immediately behind 
the reflected shock assumed to be in chemical equilibrium. 
There exist several choices for the chemical kinetic rate mech- 
anisms for high-temperature air" that give variations of ap- 
proximately ± 20V$ in the computed atomic oxygen mass frac- 
tion. The mechanism in Ref. 12 has been widely used for 
computing shock tunnel nozzle flows and is adopted here. 

Figure 3 shows computed atomic oxygen and nitric oxide 
mass fractions (the only significant air contaminants at these 
conditions) as a function of lime after flow arrival in the test 
section located 14 m from the diaphragm at the exit of the 
acceleration tube. The time ordinate in ail cases results from 



Fig. 3 Atomic oxygea (0) iod iltrk oxide (NO) mass fractions In lb« 
lest flow at Mach 17 HYPULSE condition. 



Fig. 4 Entropy -lac rtla corrrtatloa for predicting nooequOIbrium air 
compositions In an expansion tube, t* * 


conserving mass between a given particle and the diaphragm 
when going from the postincident shock state, 2, through 
the reflected shock and expansion, to the facility exit, state 3 
(Fig. I). 

Considering the holding time model results, the first parti- 
cles to exit the facility (those initially adjacent to the di- 
aphragm) experience rapid expansion rates and freeze at a 
composition near the pre-expansion value. Subsequent parti- 
cles see lesser expansion rates and recombine further, in agree- 
ment with the results of Ref. 5. For the numerical solution of 
the inertia model, five particle trajectories were analyzed. The 
presence of the diaphragm mass reduces the expansion rate for 
the first particles yielding greater recombination than the 
holding time model. Subsequent particles are processed by a 
weaker reflected shock altering the composition as shown. 
Using the pressure history from the analytical solution for the 
particle adjacent to the diaphragm, good agreement is found 
with the numerical solution there. Accuracy of the analytical 
solution diminishes for subsequent particles and they are not 
analyzed for recombination. 

Entropy Correlation 

For steady expansion nozzles the frozen enthalpy and com- 
position of nonequilibrium air at the nozzle exit correlate well 
with the nozzle reservoir entropy. 13 ' 14 Nonequilibrium compo- 
sitions resulting from unsteady expansion can be similarly 
correlated using the postreflected shock entropy, provided 
that the initial portion of the expansion occurs in equilibrium 
and at constant entropy. For this situation the expansion wiQ 
pass through the lower pressure starting conditions of other 
expansions having the same entropy. If these expansions sub- 
sequently have coincident pressure-time histories, they will 
yield the same final composition. From Eq. (4) it can be seen 
that coincident expansions have equal values of the modified 
time scale 


tm *= * iv2y (Q 

All expansions with the same r* and initial entropy will yield 
the same chemical composition. 

Figure 4 illustrates this correlation for atomic oxygen and 
nitric oxide mass fractions, ft was constructed for a range of 
initial pressures of 1.5, 15.0, and ISO MPa and temperatures 
from 3500 to 6000 K. These ranges should span the majority 
of expansion tube operating conditions for simulating subor- 
bital velocities. The correlation is accurate to within (be uncer- 
tainty due to the high-temperature air kinetics already noted. 
It applies to the particle initially adjacent to the diaphragm, 
and thus from Fig. 3, represents an upper bound on atomic 
oxygen and a reasonable estimate of nitric oxide content. 
Departure from the correlation occurs below 3500 K where the 
equilibrium starting conditions fall below the freezing condi- 
tion of expansions from higher pressure and the same entropy. 
The thermodynamic data available in the kinetics code 10 set 
the 6000 K maximum temperature limit. 
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nnir measurements on blunted cones and a 
SCRAMJCT vSSS IN HYPERVELOCITY FLOW 
I M Porter* , D. J. Mee’ and A. Pauli" 

The University of Queensland, St. Lucis, Queensland 4072. Australia 


Absind j 

Thi$ paper reports some applications of the « 
deconvolution force balance ^ , 

Simmon*. 1991) for measuring drag S , 
impulse facilities where test tunes are of b order or 
Hi*. Two basic model geometries art considered. 
Tiw fl«t ie a 5* semi-vertea angle cone with varying 
degreeso/ nose Wuntness. TUinfluena ofbtunuiess 

fn g Sion this slenderconicalMyin.^ch5.2 

air stream at a stagnation cnthalpyofH WWj « 
presented. The second application is to the drtg 
measurement on the more complex 
scramjet powered vehicle. The drag on A 
baa bon measured at sugnation enthalpies in the 

range of 3 • 14 MJ/kg. 

JJniroductiQQ 

A scramjet operating at Mach numbers above 4 is 
being proposed as the propulsion system for a 
disposable vehicle to launch small payl°^ 
catih orbit (Stalker. 1992). The forebody of the 
scramjet will be conical in geometry. • 

sharp nose would reduce the aerodynamic drag on 
ibevehicle. However, a pointed slender nose Is 
difficult to cool and may restrict the paylo . 
capacity. These practical limitations have, motivated 
research into the aerodynamic characteristics of a 
blunted cone in hypersonic (Mach number grea 
than 5) and hypcrvelocity (flight speed greater than 5 

ferns' 1 ) flow. 

To simulate these flight conditions experiments 
were performed in the T4 free pislonr JjverAo* 
tunnel facility (Sulker and Morgan. 1988). The test 
times experienced in the T4 facility art of 


• Graduate Student, Department of Mechanical 
Engineering. 

« QED Research Fellow, Department of Mechanical 
Engineering. 

« Research Fellow, Department of Mechanical 
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of 1 ms. While the drag force is a fundamental 
parameter in the design of any flight veh.e e .u 
measurement in hypersonic impulse facilities is 
inhibited by the short test time. Conventional force 
balance or accelerometer balance techniques would 
restrict the maximum model size to 50 mm in a 
facility having a test time of 1 ms (Sanderson ana 
Simmons, 1991). However, a tedmique for dmg 
measurement in hypersonic impulse »j* 

deconvolution force balance, has b«n developed by 
Sanderson and Simmons 0991). This 
takes into account the distributed mass effects 
depends on the interpretation of the transient stress 
waves propagating within the model and its 
supporting suucluit. 

The present experimental programme is 
which the technique has been used to obtain data on 
the drag of aerodynamic bodies of current interest 
One primary aim was to obtain a P rel ‘™" 1 '][ 
measurement of the influence of nose Wuntness oo 
5Sg on a slender cone in hypcrvelocity flow. 
Another was to extend the technique to measuring the 
drag on the more complex geometry of a scramjet 
vehicle Having achieved these objectives it is 
envisaged that .more comprehensive Uivesug.uon 
into nose Wuntness effects can be conducted. The 
result of applying this knowledge to the optimal 
design of a scramjet vehicle can then be measured 
directly using the deconvolution force balance to 
measure the drag in hypersonic and hypcrvelocity 
test flows. 


2 Experiments 

2 1 Facility 

The experiments were performed “» die T4 free 
piston driver shock tunnel facility (Stalker and 
Morgan, 1988). A contoured axtsymmetnc Mach 3 
nozzle was used to expand the gas from the 
stagnation region to the appropriate test conditions. 
The nozzle exit plane was 265 mm in diameter ’ and 
the nozzle throat diameter was 25 mm. The tunnel 
was operated in a tailored mode so that the static 
pressure and enthalpy would be constant throughout 
the flow lest time. 

12 Dray Measurement T tchniQUC 

The model is attached to a “sting* in the form of a 
slender elistic bar (Figure 1). The sting is suspended 


9 




Figure 1. Drag balance configuration. 


by vertical threads to allow free movement m the 
axial direction. Strain gauges located on the sting 
record the passage of stress waves resulting from the 
impulsively applied drag force as they are 
transmitted from the model into the sting. 

Following Sanderson and Simmons (1991) the 
dynamic behaviour of the model/sting combination is 
modelled as a time-invariant, causal, linear system 
described by the convolution integral. 


y(t ) = | g(t - t) u(t) dt (1) 

0 

where u(t) is the single input to the system, y(t) is the 
resulting output and g(t) is the unit impulse response 
function. Knowing the response of the system to a 
unit impulsive force it is possible to determine the 
response of the system to excitation by any arbitrary 
force via Equation 1. Alternatively, and what is 
done here, y(t) is obtained from the strain gauge 
output and a numerical deconvolution process is 
performed to obtain u(t), the time-history of drag 
applied to the model (Sanderson and Simmona, 
1991).. 

It has been demonstrated that the technique is 
insensitive to small changes in the impulse response 
function (Simmons et ai (1992)). Thus the unit 
impulse response function may be approximated 
numerically using dynamic finite element modelling. 

2 t Variable Nose BluntnCM Skndcf CfflC Mo&l 

Experiments were performed on a 5* semi-vertex 
angle aluminium cone of 100 mm base diameter* A 
total of 11 variable nose tips was used ranging tn 
nose radius from 0.2 mm to 18.0 mm in steps of 1.8 


mm. These correspond to bluntness ratios of 0.004 
to 0.36 where the bluntness ratio is the ratio of the 
nose radius to the cone base radius. The cone is 
571.3 mm long with the sharpest tip attached. The 
base of the cone was made hollow in order to reduce 
the mechanical time constant of the system (Porter et 
al, 1992). This reduced the mass of the model by 
almost 50 %. 

The 2.5 m long sting was constructed from brass 
tubing of 34.9 mm outside diameter and 1 .6 mm wall 
thickness, giving a high bending stiffness. The sting 
length was chosen to delay the interference of the 
stress wave reflected from the free end of the sting 
with the flow test time. As the speed of propagation 
of stress waves in brass is 3.6 kins* I , a sting length 
of 2.5 m allowed 1.3 ms before the reflected waves 
arrived at the strain gauges. (This estimate of test 
time takes into account the positioning of the strain 
gauges 300 mm from the model/sting interface 
(Figure 1)) 

Figure 2 shows a diagram of the model/sting 
arrangement with the finite element mesh used to 
calculate the unit impulse response function 
superimposed. As the drag measurement technique 
used here is insensitive to small changes in the 
impulse response (Simmons et al. 1992) and as the 
unit impulse response function is in turn relatively 
insensitive to the loading distribution on the model, 
the impulse response determined for a distributed 
load on the sharp cone model has been used in the 
numerical deconvolution to obtain the drag for all 
the blunted cone models tested. 

? 4 Scranye! Vehicle Model 

The model of the scramjet vehicle was made of 
aluminium and consisted of a centrebody surrounded 
by 6 evenly distributed combustion chambers (Figure 
3). The centrebody was composed of a 9* semi- 
vertex angle sharp conical forebody of length 99 
mm; a 9* intake surface extending for a further 60 
mm; a combustion duct length of 60 mm; and a 15 
thrust surface of length 55 mm. The scramjet 
centrebody thus had a total length of 274 nun. The 
combustion duct was also made of aluminium and 
had a total length of 145 mm. The stainless steel 
cowl surface angle at the thrust nozzle exit was 7 . 



Figure 2. Finite element model of sharp cone model and sting. 
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Figure 3. Scramjet powered vehicle. 


highly temperature sensitive. This made it necessary 
to use a modified Wheatstone bridge circuit to 
compensate the gauge output signals against 
undesirable temperature effects. Four strain gauges 
were used. Two gauges were altached to a separate 
piece of the sting material and placed near the sting 
so that they saw the same thermal environment but 
no mechanical strain. The two strain measuruig 
gauges were arranged so that the circuit was sensitive 
only to axial loads. 

Measurements of the pressure in the base area of 
the cone were obtained. Additional instrumentation 
consisted of measurement of the shock speeds tn the 
shock tube and stagnation pressures. These are used 
to determine the conditions of the test gas. 


The scramjet model had 6 combustion ducts evenly 
distributed about the centrebody. TJc diamc ^ 
the centrebody for the length of the wmbusfion 
rh ambers was 51 mm. The internal diameter of the 
combustion duct cowl was 69 mm and the external 
STnre S; was 71 mm. Each of the combustion 
chambers was approximately 12 mm tn radial l wi 
No provision was made for injecting fuel in 
scramjet ducts on this model 

For these tests the scramjet model wasatuchedto 
the nose of the cone in place of one of the blunt tips. 
The cone was then shielded from the flow to preve 
it from experiencing any aerodynamic forces. 

The unit impulse response function was determined 
numerically using a dynamic finite element program 
for the case of a point load applied axially at the p 
of the model. The finite element mesh for the 
scramjet model is shown in Figure 4 . The supports 
attaching the cenirebody to the cowl were 
approximate the mass of the walls of l 
combustion chambers. This was done so that the 
scramjet could be modelled with an axisymmetnc 
finite element mesh to simplify the numerical 
analysis. 


3 Results and Di scussion 

3.1 Test Flow Conditions 

The conditions in the test section were numerically 
determined using ESTC (McIntosh 1968) and 
NENZF (Lordi et al, 1966). The shock speed in the 
shock tube and the stagnation pressure were 
measured and used as inputs to ESTC to determine 
the temperature of the test gas in the stagnation 
region after shock reflection. The test gas undergoes 
a steady expansion from the stagnation region to the 
test flow properties at the exit plane of the nozzle. 
NENZF is a one-dimensional non-equilibrium code 
which is then used to predict the properties of the test 
gas at the nozzle exit plane given the stagnation 
pressure and temperature. 

Measurements made of the test section Pitot 
pressure and static pressure were found to agree well 
with those predicted by NENZF. The test flow 
properties were found to be repeatable to within 
10%. The four conditions al which tests have been 
performed are listed in Table 1. All experiments 
were conducted in a lest gas of air. 

Rlunted Tone Drat McaSU mngflU 


Instrumen tation 


Kulite UHP-5000-060 semi-conductor strain 
gauges were used to measure the system output y(t). 
These gauges have a gauge factor of 155 but they are 


Experiments on the variable nose bluntness slender 
cone model were performed at the 14.4 MJ/kg test 
condition (Table 1). The signal output from the 
strain gauge bridge was deconvolved numerically as 
described in Section 2.2 to obtain the time history of 



Figure 4 . Finite element model of scramjet vehicle and sting. 


11 





Stagnation Enthalpy 




MachNOt 


StaUa Prooauro 
hPa 


Pitot Praoauro 

Si 


Itartto Tonporvtwro 
tt _ 


Ptov Voloclty 
fcm/i 


Donaltyj 

fcfl/tnS 

0.021 


14 4 
tOJ 

jt 


SI 

5.4 

Ji- 

Ji 


It 

It 

to 

1L 


JJL 

J1L 

070 

m 


Itto 

HOP 

_31i_ 

ro. 


4.S 

4.1 

Ji 

*A_ 


0.030 

0.0»4 

im. 


Table 1. Flow conditions in a test gas of air. 


the drag on the model. Thi* drag measurement 
technique is inherently noisy ts the deconvolution 
process tends to amplify any notse present in the 
original output signal. y(t). It was l^s necessary to 
pan all the drag measurements resulting from the 
numerical deconvolution process through a 2 kHz, 0 
pole Butterworth low-pass digital filter. 


Figure 5 shows the measured drag in its filtered 
form for a nose radius of 10.8 mm in comparison 
with its corresponding strain gauge output signal 
before deconvolution. The steady flow test time may 
be identified by the regions of steady drag, occunng 
around 4.8 ms in Figure 5. 


The results from this preliminary study into the 
effects nose blunmeu haa on cone drag are 
summarised in Figure 6. At the smaller nose 
bluntnesses the effect on the total drag is relatively 
minor. The drag shows an increase from the sharp 
nose value of about 20 % at a nose radius of 7.2 mm 
(bluntness ratio of 0.144). However, beyond this 
redius the drag Increases more rapidly so that at * 
nose radius of 18.0 mm (bluntness ratio of 0.36) the 
value of drag is about 145% greater than the drag on 
the sharp cone* 


A theoretical prediction of the total drag on • 
sharp cone has been made and is also shown in 
Figure 6. The pressure drag was predicted using 
Taylor-Maccoll theory (Taylor and Maccoll, 1932) 
and found to be 159 N for the case being studied 
here. An estimate for the skin friction drag on the 
cone was made based on laminar boundary layer 
theory (White, 1974) and found to be 36 N. Tnu 
gives a total drag of 195 N for a sharp, 5 semi- 
vertex angle and 571.5 mm long cone travelling in 
air at the 14 MJAg condition (Table I). The value of 
drag measured for the sharp cone was 200 N. thus 
reinforcing the accuracy of the drag measurement 
technique used. 


^ t Drag Force on Scrami tt Vctliclfi 

This model provides a good test of the utility of the 
deconvolution force balance given the complexity 
and size of the model. AU previous studies using the 
deconvolution balance have been on simpler conical 
geometries (Sanderson and Simmons, 1991; Porter et 



Figure S. Comparison between signal from strain 
gauge bridge before (raw signal) and after (drag 
force) deconvolution for a nose radius of 10.8 nun. 



Figure 6. Drag force versus nose radius. 


al, 1992). The present model deviates from 
axisymmelry in the arrangement of the scramjet 
engine modules around the central body of the 
vehicle. All net aerodynamic forces should be axial 
for these zero incidence tests. These tests were 
performed to measure the aerodynamic drag 
associated with the scramjet vehicle design, thus there 
was no injection of fuel into the scramjet ducts. 

The first step in testing the applicability of the 
deconvolution force balance to drag measurement on 
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(*) 



Fieure 7. Comparison between (a) numerically 
oredicted and (b) experimentally measured step 
response for a point load applied to the tip of the 
scramjet vehicle. 


the complicated scramjet model involved comparison 
rfSe element computation of the steprespoose 
of the model with an experimental result The finite 
element computation was P e J fo . rmc , d v for 
impulsively applied point load at the tip of model. 
The experimental step response was obtained by 
cutting a wire attaching a mass to the tipof the m 
as in Sanderson and Simmons (1991). This simulates 
a point load applied impulsively to the tip of the 
model. Six experimental results were averaged m 
order to reduce the level of random noise on the 
tienals. The results for the computation and 
experiment are shown in Figure 7. In these plots the 
stresses are normalized with the stress that ■ 
obtained for a static application of the load to the up. 


Figure 8. Comparison between the drag time 
histories on the scramjet model in a flow of 
nominally 3 MI/k( and two static pressures. Toe 
solid line represents the drag at a static pressure of 
10 kPa and the dashed line represents the drag at a 
static pressure of 3.5 kPa. 


The results show that the finite element 
computation predicts the experiments quite well. 
The overall rise in the level of stress in the sting is 
well captured although there are some differences in 
the finer derail. Based on this result the computed 
step response has been used to determine the impulse 
response for use in deconvolution of the measured 
results. 


Figure 8 shows dmc-histories of the deconvolved 
drag on the scramjet vehicle at the low enthalpy 
conditions of nominally 3 MJ/kg of Table 1. Results 
art shown for the two different pressure levels at this 
enthalpy. The results show the characteristic peak tn 
drag at the start of the lest flow which is atmbuted to 
the starting up process of the nozzle, followed by a 
period of relatively steady overall force. This is also 
evident in the cone drag data. The period of steady 
drag force corresponds to the steady flow test time. 
Thus, it would appear that the deconvolution balance 
is suitable for testing complex models such as the 
vehicle. 



Table 2. Drag coefficients of scramjet vehicle at 
the 4 conditions of Table 1. 
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6 References 


The drat coefficients (based on the external cowl 
radius) corresponding to the four test conditions of 
Table I ranged from 0.20 to 0.33 (Table 2). These 
relatively high values are attributed to a lack or 
pressure recovery downstream of the combustion 
ducts. 


4 Conclusion 

The drag measurements made on the 5* semi- 
vertex angle, variable nose bluntncss cone model 
reveal a steadily increasing effect of nose bluntnesi 
on drag on the cone. The effect at the smaller nose 
bluntnesses is relatively small, with a 20% increase in 
drag at a hose radius of 7.2 mm. This is encouraging 
for the design of a hypersonic space plane or a 
ccntrebody for an axisymmetric scramjet where a 
slightly blunted nose is required to reduce stagnation 
point heating. Beyond this nose radius the drag 
increases more rapidly with nose bluntness. 


The results for the scramjet powered vehicle 
indicate that the deconvolution force balance can 
applied successfully to measure the drag on models 
of quite complex geometry. It is envisaged that a 
more comprehensive investigation into possible real 
effects on the nose bluntncss effects on cone drag 
will then be applied to the optimal design of the 
forcbody of a scramjet vehicle. The testing of 
various designs can then be carried out using the 
deconvolution force balance* 
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T musT measurements of a complete axisvm^tric scramjct ^ f 
IN AN IMPULSE FACILITY 
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™, sts. i; 

th« hypersonic Upu I scramjet ^onfigurmtioft. In 

integrated * xl * y,1 ^[ r force on the scram jet vehicle 
these tests the n ^convolution force balance. 

measured uslng a^ , U application to a 

This measurement technique , r< dlscuS *ed. 

complex model such at jcramjet's aerodynamic 

Result* ar. scram jet when fuel 1. 

dr*£ end the net force It 1* shown 

Injected Into the co “^* 1 hydrogen-silane fuel 

]trptflqJCTlQW 


iimraxlutely 1 as However, this Is Insufficient 
for testing scramjets. An Increase In decor "'°' u ‘ 
time Is required because valves ar* opened and fuel 
Is injected prior to the start of the test flow. To 
obtain an accurate lime history of the force on the 
scramjet these ever.ts musl also be deconvolved. The 
development of this balance and technl^i ^he 
. . „Kdbr» time-history measurements can be made 
fwlr . Tr'od “ excess" of 10 on . complicated 
Ldel has been crucial to the measurement of the net 
thrust produced by a scraajet. 

The ability to measure the overall P* rf ?™* nc * ® f * 
fully integrated .cra.jet 1. on. which should 
^ L. ih* deslxa efforts towards the develooment 
tT of e«ln*. The result, presented here 
of lh ‘* YlTst r record of the thrust produced by a 
f ully ^TntVgrVted* scraajet which ha, been tested In 
an Impulse facility- 


Extensive 1 nve s 1 1 8 *^ 1 “"^^^ted to* determine If 
experimental, *»ve ' thrust to be 

: C c^c S .Uy P vUbU -^rsonlc night . ^hough 

scramjet 1. "'tt l, to say. flight 

Integrated eonflguraU • lnUkes> combustion 

tests of a scramjet exterior surfaces, 

chambers, thrust s ^ f4C * lw0 way5 . A scramjet 

-SSTaimS;^ mum;-- - *• 

- tr-.vys ? 

test the ,X,S *“T? t produced by the scramjet was 
The drag or net l ^ ust d / co „ vo i u uon force balance 
obtained using i toot 11 The time-history of 

(Sanderson and Simmons Cl ’ , obtained by 

the force acting on the „ dt of the 

numerically deconvo v ' ^ t0 the scramjet. 

mxlel strain In a «tl*»» or i.inally designed so that 

rarsffTssss /... . m— - 


ICOEL DESCRIPTION 

A schematic of the scraajet. fuel tank and * t ! n *,!* 
A *"', n f «„ r , i. The scraajet Is connected to the 

fl£l tank which la tern Is connected to the *tjh»- 
The scramjet Is -cylindrical In shape with a conical 
J^e^v and thrust surface. It has an mls-lrrti- 
t TbSdy and a stainless steel cowl. The fuel 

tank Is ^cylindrical In shape and Is 
1 , . , /...i tv. s tlng Is a brass cylindrical 

: « ...... ... 

testing approxiulely 200 mm dewnstre.a of the 

The'st lng and model are suspended by fine threads so 
thlt the combination Is free lo move In the 
direction of the axis of the tunnel. The fuel lank 
and sting are shielded from the lest flow. This 
is required lo ensure that the forces 
which a^e measured ar. only the a.rodyna.lc forces 

The' over all** length of the scramjet Is 274 mm and Its 
external diameter Is 71 am. The conical forebody 
h« ^ degree half angle and the thrust surface ha. 
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around the cyH^Ual cenUebod/^ 1* -.busUon 
fiction. Se-trebody * cylindrical -Uh 
dilt eler of 51 ** thgfc ^*1 and the fuel 
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tank is 75 »« 1°"*- t , c 55 •* downstream of the Vhe 
leal* on this P'P* *1*^ \erodyna.lc drag on this de 

s-* s«rs^ p- - - d - on ihe a 

ruel" could be Injected up«lr.» ■ °\£J 0 °T ihe «j 

fontr.t }S3i ^tr^' 'Al- - 

J 3 S 5 f ^ted^t 1 SI - 

2ES t°o n ^ -4 tlon c^* "* upslre " - 

injector* were plugged . is one injector. P r 

In each co.bust Ion <ha.l*r there M degree « 

Each injector U » ! Z “ scrwjet and 

to the centreline © co.bostion chamber. < 

centrally at the froa th e fuel tar* J 

£*EV pr.»- * 

S5i.“f “f fs'^r.rf. ssrts s 

MSI flow rate of the fuel- ^ » as s flow ( 

flow through the injec changing fillip 

rate of the fuel ls * trie ^ 

pressure of the fuel la t fast acting valve 

Fuel flow 1* controlled F fy^t tank. The 

located between the *craaje {depending 

JST«.»» *•“.*, \yz 100 

on filling pressure) ana c 

»*. 

•pnr mure HIASUPDOCL jigg^ 

The force balance .us'd Ko Sanderson and 

thrust 1. describe In graining the 

Simoon* (1991) » acting on the test aodel 

Uae-hlstory of the force |l^g ^ # long sllng 
bv nonltorlng the »**J l aodel. When the now 
attached to the tase of «* % omn , aerodyna.lc 
starts around the t0 propagate and reflect 

forces cause * lr * s * . l>nk . These waves are 

within the aodel •«* f “*' tnd , re altered by the 
transmitted Into th ”* #|> C an be aodelled as a 

strain gauges. * y ^l by an Input u(t) - the 

linear syste. characlenaed «y y(l) - the 

axial force on the »od i .pulse response 

strain - e »* ure * *\ l JT can be written using a 

g (t) relating the two. inis 


convolution integral. 


y(t) - fg(t-r) u(t) dx. 
’'o 


— « s’n-Lsrr r^x\”,r£ 

axial for '* ^esswnsT^whlch character lies the system 
The ‘“P uls * J? ^. t her experimental ly or by using 
*•" ^ Hnlte element analysis. For the present 

dynaalc ully determined l.pulse response 

tests an expe y found by suspending the 

h *J ,^",1 bTi wire attached to the Up 

«odei »nd '^8 * .. f lh * # cramjel. The wire 

of the conic* l ore J of the forebody to 

„„ then cut close to # P^,. loid This 1 . 

produce « sudden re j om( j 4pp ned at the 

equivalent tom * € ^‘ resulting strain gauge output 
tip of the »odel- The res i ** the #t ep 

(figure *» ~” ur * ‘Jl ‘V “W* response I. 

^r?ound by differentiating the step response - 1 th 

respect to tlae. Billiseconds 

The opening of the f«-l «» ^ flotf , d ded soae 
before tbe arrlv ^ fore€ balance seasureaent 
coapl ca on . v | ous jy encountered. Previously 

USr£n and ll^ns ( 1991 )) strain gauge signal. 
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r. Strain gauge output af a function of III 

fir "unit step load applied to the scraajet nose. 
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had onlv been processed to the point In tl»e * 
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ijr^w-ss narssL* -2 « 
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^deconvolved through t number of stress 
reflections. 

run. description 


Initially, hydrogen was used as .lule^thrMt **«•« 
« I'** ^is W«J th.‘t U t e he co-bust.on 

chambers are too short for the fuel to fully coabust 
al the conditions targeted. To reduc * 

combust Ion length, silane (S1H 4 ) was added to the 
hydrogen. 13X of the fuel's volu * e “ a * *‘i 
Fuel was Injected at an equivalence ratio of 0.8. 

F l t^!cal record of the fuel pressure Is given n 
fl u r It can be seen that the fuel pressure Is 

constant during the test tl.e. ' 

reialr.s constant for approximately 3 as. after ^ h 

falls rapidly as the fuel Is drained from the 
r .i .ink It is critical that the fuel valve is 
opened so that this 3 as period Is coincident with 
Vest flow. Ideally, the fuel pressure should 
J«ch It. saxlr.ua approxl.ately 0 5 as before the 
II - 1 flow arrives. If fuel Is Injected early then 
the Bass flow rate of the fuel will not be steady 
and furtheraore. the tlae-hlstory of the strain 
induced into the sting fro. the Injection of fuel 
■av not be fully recorded or will be excessl ejy 
,<4. Consequently, the force on the sera. Jet coul 
not be properly deconvolved. 

RESULTS 



Figure 3. Fuel pressure as a function of tine. 



Time (ms) 


Mg r S the 4 - tesU 6 ‘VST X TTrn^el 

represents* the VImms t.Arlng the tests. The flow 
arrives at the test section at t ■ 0.5 as. 

FLOW PROPER T ) CS 

In the operation of a reflected shock tunnel the 
list gas H first shock heated and Is then expanded 
through a contoured -rle^The U.pe^.ture of the 

expanded* l/^deteralned using ESTC (McIntosh 0 968)1 
to be 2700 K. Properties of the test g»s l 
lest section are then determined using HOIZr ( Lordl 
et al. (1966)). In the experiments described here It 
I, assume that the test gas Is In ch.alc 
equilibria as It expands down the nozzle. It has 
Sen deter.lned that. In the test section the test 
Hs has a velocity of 2400 a/s. » density « X I «• » 
kg/a . and a temperature of 350 The 

enthalpy of the flow 1, 3J HJ/V, » “ 
equivalent to a flight speed of 2600 a/s. 

F leu re 4 displays measurements made of the static 
and Pitot pressures. It can be seen that e ** 
flow establishes Itself approximately 0.5 « »f er 
the arrival of the first shock. Sub^u'ntly. the 
flow is quasi-steady for approximately 1.5 as. 


Figure 4. WaJJ end Pilot pressure ms m function of 
tine. 

FORCE MEASUREMENTS 


Figure 5 displays the ou'put fro* the strain gauge 
mounted on the sting. The thrust produced by the 
Injection of the fuel Is observed as a relatively 
slow oscillation prior to the arrival of the test 
gas. When the test flow arrives it can be observed 
that the output is more complex. 

The output In Figure 5 Is deconvolved using the unit 
Impulse response function to obtain the time-history 
of the net force on the scraajet. The result of 
this deconvolution is displayed In Figure 6. 

Figure 6 displays three such deconvolved 
tlme-hlstorles. These are the net forces on the 
scramjet when 

(i) fuel is Injected Into a test gas of air, 

(II) fuel Is Injected Into a test gas of nitrogen 
and 

(III) fuel is not injected and the test gas is air. 
In the latter case It can be seen that the scramjet 
experiences a drag-history which basically reflects 
the slight decay observed in the Pitot pressure 
(figure 4). It can be seen that the drag Is 
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approximately 140*25 K. 

When nitrogen is used as a test gas and fuel la 
Injected the drag on the scramjet Is approximately 
90*20 N. Hence, the drag Is approximately SO N less 
than the fuel off test In air. There Is a slight 
difference in flow conditions between the air and 
nitrogen test gases, however, the drag Is reduced 
primarily due to the presence of the fuel. When 
fuel is injected. It can be seen, prior to the 
arrival of the test flow, that the fuel Induces a 
thrust of approximately 30 N. It Is speculated that 
the the remaining 20 N reduction in drag results 
from the different test gases and fuel air 
interactions which possibly change the skin friction 
coefficients within the combustion chambers. 

Figure 5 also displays the drag and thrust on the 
scraajet when fuel is Injected into air. It can be 
seen, as above, that prior to the test flow the fuel 
Induces a net thrust of approximately 30 N. As the 
flow establishes itself across the model the drag 
Increases, following closely the result obtained 
with the nitrogen test gas. However, once the flow 
Is established Ignition of the fuel occurs and the 
drag Is substantially reduced. 200 ps after flow 
establishment the thrust created by the combustion 
of the fuel Is sufficient to balance the drag. I os 
after flow establishment and during the test time 
the net thrust is approximately 50 N. During the 
test time the net thrust peeks at approximately 95 N 
and is always greater than 40 N. 

CONCLUSIONS 


The Impulsive force balance and the deconvolution 
techniques described in this paper are sufficiently 
advanced to measure the net force on a complicated 
model . 

The use of an experimentally determined unit impulse 
response function has made it possible to determine 
with good frequency response the net force 
time-history on a model such as the scramjet, where 
the time-history of the forces applied to the model 
are of the order of 10 ms. 

The scramjet used In these tests, although producing 
more thrust than drag. will need to be more 
efficient to become viable. The technique described 
above has made it possible to test new designs In 
hypersonic Impulse facilities. Thus, definitive 
values for overall performance of scramjet designs 
can now be obtained relatively simply and cheaply at 
hypersonic conditions. 



-3 -t -1 0 1 C 

Time (ms) 


Figure 5 . Strsin gauge oufput as a function of 
tine. Fuel Injected into sir 



Time (ms) 


Figure 6 . Net force on the sersejet ss s function 
of tine when fuel Is injected Into sir and nitrogen 
test gsses and when fuel Is not Injected snd the 
test gss Is sir. 
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introductio n 

This note reports tests in a shock tunnel in which a fully integrated scramje, configuration 
produced net Uuust. The experiments no, on,y showed Urn, imp„.se facilities can be used 
for assessing duns, performance, bu, a, so were a demonstration of the app.icton of a 
new technique « to the measurement of Onus, on scan*, confutations m shcck 
tunnels. These two developments are of signiftcance because scramjets are expected to 
operate a, speeds we., in excess of 2 km/sc. and shock tunnels offer a means of 
generating high Mach number flows a. such speeds. 


THF, MO nFI AND TEST FACILITY 

A sketch of the scramje, model and fuel tank is shown in Figure 1. The scramje. 
centrebody. shown in streamwise section in the ftgure. consisted of a conical fotebody 
with 9- half angle, a cylindrical section of 51 mm diameter, and an afterbody of .0* half 
angle. It was party surrounded b, a axisymmetric cowl, which had an internal diameter 
of 67 mm over the parallel section of the cenoebody, and was of 71 mm outside 
diameter. Filler pieces, which are no, shown, divided this parallel section into six 
constant area combustion chambers, each of which subtended an angle of 26* a, the 
centreline. Fuel was injected through six oriftces, each 2 mm in diameter and angled a. 
30- to me centreline. They were located a, the upstream end of the combustion 
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chambers. The filler pieces between .he combustion chambers extended ups«am in ttre 
form of intake compression ramps which processed the flow in .he forebody shock layer 
through two shocks, each of 8' deflection. The leading edges of the cow. were shaped to 
prevent these shocks spilling from the intake, thus forming a convergent duct leading to 
the combustion chambers. The radial dimension of this due. was 10 mm. 

The fire, supply system is also shown in the figure. The fite, tank was filled before a test, 
and the supply valve remained closed until the test was initiated. The recoil of dr. shock 
tunnel closed a switch to operate the solenoid valve, which opened me supply valve. As 
the fuel flowed to the model and fuel injection orifices, its pressure was monitored by a 
PC B piezoelectric pressure transducer. The system was designed so that sonic flow 
occurred a, the orifices. As shown, the fite, tank and valve assembly was shielded from 
the flow by an aerodynamic shroud. In normal operation, the upstream end of this 
Shroud was located 50 mm downstream of the termination of me scramje. afterbody, and 
,o check mat mis was sufficient for me shroud to have no influence on me forces on me 
scramje., a test was run with an insert to increase mis disurnce to .25 mm. No change in 

the measured force was observed. 

The experiments were done in me free piston shock tunnel T4 a. me University of 
Queensland which had a shock tube 10 m long and 76 mm in diameter. A mixture of 
Argon with Helium was used as the driver gas. and me rario of me rwo was varied as me 
stagnation enthalpy was varied to ensure that a sufficient period of approximately constant 
nozzle stagnation pressure was maintained throughout me tests. The contoured 
hypersonic nozzle used for me tests had a throat diameter of 25 mm and a tes, section 
diameter of 250 mm, and was operated with a reservoir pressure of 37.5 ± 2.5 MPa. 
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THE STRESS WAV E FORCE BAL ANCE 

A stress wave force balance, as original outlined in ref . by Sanderson and Simmons, 
was used ,0 measure the aria, force on the model. They employed dm melhod ,0 
assure dm drag on a shun cone, bu, 1. bas since been develop* for longer cones 
and for long cones wUh me load distribution imparted b, a blun, nose ». Essentially, « 
involves measurement by a strain gauge of the time variation of stress induced in a stress 
wave bar by me forces on me model, and deconvolution of me resultant record to obtain 
ne, force on the model. As shown in Figure 1. me stress wave bar was attached to 
me downstream end of the fuel tank, and me strain gauge was mounted some 250 mm 
from the junction. The stress bar was free a, me downstream end and was suspended by 
wires, which did not influence the propagation of stress waves. 

The relation between the force input u(t) and the strain gauge output y(t) can be written 
by using the impulse response g(t) relating the two. 


as 


y(t) = J o ‘ g(t " t) u(t) dr 

Dynamic finite element computer simulations showed ft* me sum of the forces 
distributed along the length of me mode, could be accurately represented by a single force 
applied a. the nose, so me required g(.) was obtained experimental,, b, venically 

suspending me mode, a* strcss wave bar by a fin. wire from me tip of me forebody, 

, no the wire The resulting strain gauge response was then 

and suddenly severing the wire. 

differentiated with resfWCt to time to find g(t>. h - — ■» to establish g(„ for a 
period of several millisecomfs, in order to accomodate me generation of thnts. by fuel 
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injection before the initiation of the test flow. As the stress wave bar was only 2 m long, 
and was made of brass, this implied that the impulse response included stress wave 

reflections in the bar. 

A computer program was prepared to deeonvolute the above integral with this g(.) and 
applied to yield the axial force for each test. 


F.YPER^fF.NTS ANP RESULIS 

The experiments were conducted with a fuel consisting of 13% silane (SiH.) and 87% 
hydrogen by volume, yielding typical results as displayed in Figure 2 «. Figure 2(a) 
displays UK pitot pressure and static pressure of Uk free stream measured during a rest, 
showing that the flow reaches an approximately quasi-steady state about 0.5 milliseconds 
after fust arrival at the test section, and that slate subsequent persists for approximately 
1 millisecond. Records of the axial force experienced by the scramjet are displayed in 
Figure 2(b). It can be seen that with no fuel injection. Uk scramje, experiences a drag 
force of 140 ± 25 N. When nitrogen is used as the test gas. and fuel is added, the drag 
on the scramjet is reduced substantially. The test conditions with nitrogen differ only 
slightly from those with air. so most of the difference is attributed to the thrust due to 
feel injection alone. When -he fuel is injected into air, it can be seen that, prior to 
arrival of -he test flow, the feel produces a thrust of approximately 40 N. and as the flow 
about UK mode! establishes itself Uk drag increases, closely following Uk nitrogen result. 
However, once air flow through the combustion chambers allows ignition of Uk fuel. Uk 
tfeust due to combustion increases until, during -he test time, it is sufficient to produce a 

net thrust on the scramjet of more than 50 N. 
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S CRAMJET THRUST MEASUREMENT IN A SHOCK TUNNEL 

TV variation of the thmst whh tunnel station enthalpy, or with the computed test 
suction velocity is shown in Figure 3. I. can V seen that ahhough the net thtus, is 
positive around sugnation enthaipies of 3 MJ/kg, it falls off rapidly with increasing 
stagnation enthalpy. It is though, .ha, this is largely due to the dissociation of combustion 
products in the combustion chamber. An attempt was made to increase the maximum 
Utres, measured by either increasing the amount of fuel injected, or by lowering the 


tunnel sugnation enthalpy. In 
onset of thermal choking. 


both cases unsteadiness in the flow resulted, suggesting the 


In order to confirm tha, the flow in the combustors was supersonic when combustion was 
taking place, pito. and static pressures were measured a, the downstream end of one of 
the ducts jus, upstream of the expansion caused by the upstream comer of the afterbody. 
At a sugnation enthalpy of 3.3 MJ/kg the ratio of Urn two was 4.1. yielding a Mach 

number of 1.75 ± 0.1. 


The performance of me scramje, configuration tested was no. impressive, as the 
maximum observed net duns, of 60 N corresponded to a fuel specific impulse of only 200 
seconds. However, no serious attempt was made to optimise the configuration for ne, 
thrust, as is witnessed by a fire, -off drag co-efficient of 0.13 a, a sugnation enthalpy of 3 
MJ/kg. Subsun, ial improvements may be expected from proper attention ,0 aerodynamic 

design. 
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FIGURE 2. TEST RECORDS SHOWING GENERATION OP NET THRUST 
(STAGNATION ENTHALPY = 3.2 KJ/KG, FUEL/AIR 
EQUIVALENCE RATIO - 0.90) 
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•n, r .„ measureme nt in a M) scramieU L onlf ■ SwnTu lt h 

/ ??iT I/??'? 7 


13-0^ 

4-93>S>0 


Supervisors: J.M Simmons, D.J. Mee 

Th* two^insensionaHhi^t^OTzle^pr^nfc^ Seles. A 

axisymmetnc as in the case accommodate the internal now and achieve some 

Iwo-sting system has been chore * 300 mm long and together the nozzle 

symmetry. ramp waits is 1 1" and the area ratio is 4.76. 

and slings weigh 6.55 kg. B ^ t he small ramp angle and the internal 

The situation is complicated by the fact that wm, ^ ^ # ^ ^ |hnls , 

pressure on die norz'e walls. P^> y aU fa the stings). Although bending 

which is to be measured (i.e. the tensi P F 6 s waveSi thc system 

stress waves travel at most t^'L ^e econ ^ng wL originally only used to preserve 
needs to be stiffened against bending. The ^ccond g & different at some 

synuneuy. However. 

ssrsri-r-- 

distorting their cross-sectional wouldno^ significandy alter the propagation of 

Finite element analysis showed that th s wou a no g ^ 

the axial stress wave in the sting, w 1 e e ^ f ^ nozz i c t0 a step point load 

„ jshs zszssi- s. - ~ — — - * 

mm Bp a. die exit of die combustor will “ s " Ler of die lip there is 

allowing free 7^"' ° cc ^'““di the nozzle. This is show^h Figure 3. The alignment 

*- ^^^"e .OP surface of the nozzle samp wads in order to keep 

them out of the nozzle exit plane (see Figure 3). of ^ nozzle and a pitot rake 

There are static pressure lappings m the ^ramp wtdBQ^ ^ ^ ^ 

provides a survey of the flow at t e exi o _ removed in order to measure the pitot 
pressure'acmss thTcornburforexi,. ^The scramje, is 600 mm long and fuel is injected from a 

“" tral Thethrust'measured via the strain 

Should be less than the thmst calcu ate jom s^ ^ pressures in the nozzle are different on 

to the skin friction. g q MJ/ke The open and closed triangles represent 

either side for the non«ombust,on ca*s at 9 MJAg. J ^ respectively. The 

the static pressures on one thrt, Me mdkaled by die open square 

corresponding pressures on th pp different pressure on the first half of each 
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. * — -ns a : sssxr sanssr- 

5SSES!*-' ~ » 2SKSS 

SS.cc:s^;-‘-“ 

farther investigated. ilBide rhe nozzle shielding have been measured 

far aST^ *‘ w» ta “ •“ re8ions °"' y ^ 

approximately 3 milliseconds after the test ; tim . ^ ^ resu )t shown in Figure 7, 

" A Mach 4 contoured nozzle, snsedm to on^df the tost 

«* ‘ hro5 ‘ calcula, j itf'Svolved strain measurement of the net thrust for the cases 

surfaces is compared with the decon approximately 9 MJ/kg nozzle supply 

zzgTJ : C “no 1 <*» — « «—> - - w “ iOT - - m 

thrust due to combustion is visible in this result. 


Thrust Nozzle Stress Wave Balance 
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Fig- I 





Connection between Scramjet Combustor 


and Thrust Nozzle 


Support Wire 
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Static Pressure (kPa) 


Static Pressure in the Thrust Nozzle 



Distance along Nozzle 


Fi&. 4 
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Raw Strain (in Volts) 


Raw Strain Gauge Output from each Sting 
Pressure Loading Asymmetric 



Time (seconds) 


Fig- 5 . 
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Deconvolved Net Thrust & Pressure Thrust 
with and without Fuel at 9 MJ/kg 



Fig - - 7 
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SCALING IN A SCRAMJET COMBUSTOR 

M.V. Pulsonetti and R-J- Stalker 9 ^ 

Studies of scramjet combustor ““'“l? '"s^ obatarf n in a "large^combustor of 

in ,he and in Fi8Ure ' lhe of fc *" 

combustors varied by a factor of 5. 

„ was anticipated Urn. the MS SM- 

length if*, temperature, the ^"^J^whils. showing for the necessary 
the equivalence ratio were the sam • ^ ^ different shock tunnel nozzles. The 

difference in pressure die two m _ g nozz i e with an intake consisting of a pair of 

large model was tested I m a "™“ , M hi * h pro du'ced shocks to reduce the Mach number 
opposing wedges preceding the t mod . J u was then expanded to 

,o P approximately 3A before the tow «UemdJhe mo , „ omM M , 4 nozzle, 

*£££ «**• — —*■ -*»■ 

The wedge intake for r the tlarge stagnalkm Sthalpy, it was possible 

Mach number for both models. Th , y ^ aream ^ ilion was matched 

also to match the temperature and > **T*^ 1 ' lhe n02 zle throat, and the initial 
by operating at the same a , ^ chemica l -freezing- took 

nozzle divergence angle were sim , goring that this was not changed by the 

&T.te“*e large midel, it was possible m turange that the model flows 
mcperienced approximately the same free sueam composmon. 

The two combustors were of constant area foT'so'x'lOO rrun.'and the small combustor 
large combustor was of la^e l^ion of each duck and was 

was 10 x 20 mm- ^ in J ^ Unfortunately, manufacturing difficulties made it 

ES& t £ ££- 7'““^ tLrrr: 

f— *• ■*«“ «* “■ differen “ reU,ive 

injector thickness should not be significant. 

Wall pressure distrib “°“ tbe P ressure r asure ‘\. at JJ 

trttx sr — ta 

with respect to the height of th • . .. DreS sure distributions, confirming the 

. - a- - — - 

length. 

The flow conditions for , t ^ expe ? me 'J^ w^T ^x^O^o^both ducts. Therefore the 
the Reynolds’ number, based on duct length, was 3.7 x 

mixing layer was turbulent 
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, „„ at lnW( . r etannation enthapies the sudden pressure rise seen in 
When le *^ rt ^"Scaling and IgniL Effects in Scramjets" by M.V. Pulsonetti Figure 
previous tests (e-g- *^. 57 ^ Supplement 8 , 1993) was in evidence in both the 

LeTd tSe small ducts. This also scaled approximately as the product of pressure with 
length. Further analysis of this data is under way. 


H, m = 7.4 MK/kg 



Large Scramjet 

Small Scramjet 

p 

18.6 KPa 

91.2 KPa 

T 

1500. K 

1500. K 

M 

4.36 

4.30 

u 

3400 m/s 

3200 m/s 

4> 

1.36 

1.27 
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Abstract 

Experimental results are given for the *P* clf ^ 

rsjc =s «*„ •« 
EE s £■» . 

and model configuration. 

Jnt rodugt lOB 


chamber duct length Is finite there is insufficient 
time for the reaction to take place. Hence, the 
specific impulse drops to zero simply because the 
fuel has not burnt. 

At higher temperatures the decrease in specific 
Impulse does not primarily result from incomplete 
combustion due to Insufficient combustor length, but 
is the result of a decrease In heat released. As 
the free stream temperature increases the percentage 
of water formed when the reaction reaches 
equilibrium, decreases. Therefore, there is less 
heat released at the higher temperatures. As a 
result, when the freestream temperature rises there 
is less energy available for thrust production and 
therefore the specific impulse decreases. 


In the past ten years there has been renewed 

Interest in air breathing engines for the purpose o 
Interest i <w of the ®ore promising 

hypersonic g * SC ramJet. Many computational 

o . *e ramjet’ » efficiency have been 
«^ (« « Ker rebrock (1992)). however, there » 

limited experimental data available which can be 
!S“ validate these predictions^ One of the 
our noses of the experiments reported here was 
SbUln good experimental data which can ta _“*"**• 
wlldatm codes which predict the overall performance 
Of a scraajel. 

lb. overall performance of the scramjet 1 . 
measured In terms of specific Impulse. It Is 
accepted that specific l-pul««. 
at 1000 . are required before scraajets can 
. Ki. alternative to rocket propulsion. 
Experiments have been performed by SUlker and 
!£^>n (1984) where the Mach nu.ber *[ y “£ 7 ^ 

the* combustion cha.ber was nominally 3.7. They 
showed that the specific Impulse peaks at a flight 
sneed of 2 2 km/s and decreases rapidly to zero If 
twilight 4 eed is less than 2.0 km/s 
mono tonl cal ly with Increasing flight l * 

zero at 4.0 km/s. Although the *^T| 

■ore relevant parameter to hypersonic fll|W. 
temperature of the test gas Is probably more 
relevant to the chemistry of combustion. In the 

exnerlments reported here the flight 4pe 
S^ctW relate to the test gas temperature. As 
the flight speed Increases then so does the test gas 
^rXur/cSe* table 1). The flight speed 
whtch the specific Impulse peaks In Stalker mn 

w '% •a^rsrsji 

approximation proposed by Pergamont (19631 1 for the 
- of a hydrogen-oxygen reaction, it 

uTiArti- P.M- «-» 

c 1>r „ r or temperatures less than 

5S21SK X‘\ iw." .. «- «— «- 

g |W W «• — * — - - — 

Im All r» «M‘ 


The scramjet is most efficient when the 
temperature at the intake to the combustor is the 
lowest which can produce combustion within the 
length of the combustor. This means that for 
reasonable size combustors the temperature at the 
Intake must be maintained near 1000 K. Therefore, 
as the flight speed increases. if the peak 
efficiency is to be maintained, the Kach number at 
the entrance to the combustion chamber must 
increase. 

The main als of the experiments reported here 
was to determine the specific impulse as a function 
of flight speed (or temperature) at a Kach number 
greater than that used by Stalker and Morgan. A 
Kach number of nominally 4.5 was chosen. 

Experiment Description 

The experiments were undertaken in the free 
piston shock tunnel T4. Flight speeds range from 2.S 
km/s to 5.5 km/s, the Kach number at the Intake of 
the combustion chamber was nominally 4.5 and intake 
static pressure ranged from 30 kPa to 120 kPa. The 
tunnel was driven so that the static pressure and 
Kach number were steady through out the test time. 
The test time ranged from 1 to 2 ms depending on the 
test conditions. Flow properties at the intake to 
the combustion chamber were monitored by measuring 
the static and Pitot pressures. The remaining flow 
properties at the Intake were determined numerically 
f roa the non-equilibrium codes KEKZF (Lord! 
et. al. 1966) and ESTC (McIntosh 1963). Table I. 
displays the conditions at the intake to the 
combustion chamber for each experiment as predicted 
using theses codes. The test gas was air. 

The scramjet consisted of a combustion chamber 
followed by a li° divergent section as shown in 
figure 1. The scramjet combustion chamber had a 
rectangular cross section of dimensions 27 mm x 54.5 
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Figure t. Schematic of scram jet combos tor and 
thrust surface. 


M The combustion chamber was 576 « long. The 

i':rr st run. Ttrrs? «« 

;L XitU .!»»■ ih , re « i ’, * 

* , tv. leading edge of the Injector - 

Of the entrance to the combustion 

chamber. 

static pressure measurements were taken 
II. c!“r.n~ of .« -“S 

%z *is. rs ;r=oiX'^ >"7« 

KXS" “ <««• ff» 5 

|ts. 


AT 


fe - * ) 


( 2 ) 


Experimental Analysis 


direction was determined fro. P It was 

measurements made on the each location 

assumed that the pressure «•«**£ 

- T» S ?,ST: «3- en dpoint at 

SSTST ~.s“e-nt - -£ S^xVo^on 

thrust va, obtained * l *? * #l £ pressure 

the thrust surface t£^Ld StiSlnlh. the 

sssr? ; ir»r riru- «— 

(l.e. sin(ll°) x total force). 

The efficiency ofthescrsajet 
terns of the specific lspulse. 

Impulse, 61. Is defined as 


41 . tT/m/t 


(1) 


„ at u the trust Increment produced by the 
here AT Is tne yu»* fc lg the mass 

arssniw-^w. 

he thrust produced *^* en f f l j, not Injected 
To. the -thrust* produced when ^ the 

[fuel-off). n* 1- * 1 fuel-on and fuel-off runs, 

renditions be t “ e *" th ^/j.tlon In conditions the 
fo take account f »v_ pitot pressure and 

thrust 1. first ®or**l l, *db^ normalized thrusts Is 
the difference t*t«e« dimensionality this 

obtained. To obUin y* . lhe p ltot pressure 

difference is then multlpU** by 
fro. the fuel-off run Thus. 


The flow properties displayed In Table l » r « V** 
flow properties of the fuel-off rtm. ow 

properties of the fuel-on run were within 10* of the 
Tuel-off properties. 

The incremental specific impulse is an estimate 
of the specific Impulse produced by a complete 
scraajet. That Is. a scramjet which has an Intake, 
in addition to a combustion chamber and thrust 
surface. This estate Is exact If there are no 
loses due to friction and the wave drag Is 
to zero. The Incremental specific Impulse Is • g°?j 
estimate If the frictional loses are small and the 
Intake, combustion chamber and thrust surface are 
designed to alnl.lze wave drag. 

Experimental Results 

Figure 2 displays a typical record_ of the 
static pressure during the duration of the flow. 
Fuel was Injected at an equivalence ratio of 
This particular record Is the pressure measured^ 
mm downstream of the corner. The test time has been 
marked.. 

Figure 3 displays a typical distribution of 
pressure during the test time .at the e o 
combustion chamber and along the t^ustsurfac.when 
fuel Is Injected. The distance displayed on the 
x-axls Is the wetted dl.tmncm downstream of to 
combustor Inlet. The beginning of 
surface occurs at S75 «. It can be seen that there 
la a drop In pressure downstream of this point 
followed by a rise and a steady decrease to zero. 
The rise Is also present In the fuel-off case, 
however. It 1. greatly reduced. The primary cause 
for this rise Is the conbustlon of the fuel. 

Figure 4 displays the thrust as • function of 
time for a typical fuel-on run. It can be _ seen that 
the thrust slowly rises then plateaus (to within 
5X). The slow rise Is reflecting the traversal of 
the shock down the thrust surface. The trust 
measurement Is taken In the plateau region. 

| a the experlnents presented the Incremental 
specific impulse was determined as a function of 
flight speed for both hydrogen and ethane fuels. 

Hydrogen Fuel Results 

The hydrogen fuel was Injected at mass flow 
rates which produced equivalence ratios of a half 
and one (stoichiometric). The specific impulse was 
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rigor, 2. Typical sf tie pressure record. 



Figure 4. Typical thrust record. 



Figure 3. Typical static pressure 
the end of the coatustlon ch saber 
thrust surface. 


distribution at 

slow the 


w* ml each flight speed for two nominal Intake 

display* the result* of these experiments. 

it can be seen froe figure 5 that for an 
. ‘L, ™ ratio of one and at the high pressure 

equivalence ratio # of the resu n* are 

condition the i g ‘ Stalker and Morgan. 

.1.1 lar to ^f'^yeTon figure 5. The aajor 
which are s\so P * ^ specific impulse has 

decreased fro. 1800 s to greater flight 

the peak ‘^ti.heXch wibTr. The free.tre.m 
t^raUri at the peak specific 1-pulse at a Hach 


of 4.5 is approximately 800 K, which is 

itthtly lower than that recorde<s ** Si *}i?L *** 
t) for the lower Hach nu^er condition. 

These results are also qu.llt.tlvely con.l^ent 

- t-sr r szr * sr- fs 

sses-a sr^sjsrssi. “ * 

At t>e lo- pressure condition It can be seen 
hat at „ equivalence ratio of one the sharp 

B the ‘^i t C io^olV 

,v - 

rlther combustion or mixing of the 
yeeur. 

At an equivalence ratio of a half the trends 

ss«t Si-iss- •» * g KiSStk 

ratio is generally higher. This Is * j 

« "n 1 . 1 t:f- » rr,’~ .?«.»••«• 

„ SSJw ...... Li. « >• • 


Ethane Fuel Results 


Figure 5 also 
produced by ethane 
pressure conditions, 
the same Injector 

experiments. 


displays the specific impulse 
fuel when Injected at the high 
Ethane was Injected through 
used for the hydrogen fueled 


The fuel vas Injected at a -ass flow rate which 
»duce an equivalence ratio of a half. It can be 
tn that In this -odel a relatively small «P«clflc 
pulse is produced by the ethane. 
evlous experiments presented by Pauli (19921 
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5 specific Inpulse *s a function of flight 
s£* *t conbust ion chLber lntkke H-ch renter of 

I'lonbustlon chtnber intake Kach timber of 3.7. 


5.5 


indicate that the length of a cortustor required for 
.thane fuel Is considerably longer than that for 
tydrogen. Hence. It la not surprising that a poor 
specific lapulse was obtained using this aodel. 

Exclusion! 

The change In specific 1 .pul ..with fillet 
■peed at a co.bu.tlon lntakeK.chnu.ber of«.5 .how 
similar trends to that observed at lower "■?* 
22S are qualitatively sl.llar to the trend, 
predicted using kinetic theory for a hydrogen fuel. 

flvdroaen fuel produces a greater specific 

!r^..’“ur..r^uT , v h rru u ^ 

Incomplete combustion of the ethane, 
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Notation 

9.8 m/s 1 
Total Enthalpy 
Specific lapulse 
Mach Number 
Hass Flow Rata 
Static Presstjrs 
Pitot Pressure 
Thrust 
Temperature 
Gas velocity 
Equivalent Flight Speed 
Equivalence Ratio 

H2 Equivalence Ratio which is nominally ooe 
B2 Equivalence Ratio which Is nominally a half 


Subscripts 

po Fuel-on Properties 

off Fuel-off Properties 
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H 


V 

Ks 

T 

km/% 

MJ/kg 

K 

2.8 

3.9 

700 

2.9 

4.3 

870 

3.3 

5.3 

1020 

3.5 

6.2 

1230 

4.2 

8.8 

1590 

4.9 

11.8 

2370 


P P 
kP» kg /* 3 

tow Pressure Runs 

35 0.189 4.5 

30 O.U 4.4 
28 0.09 4.5 

30 o,09 4.4 

35 0.07 4.3 

39 0.05 4.3 


U # :s 

km/% 


2.4 0.9 3 42 

2.6 1.0 

2.8 1.0 3.^3 

3.0 1.0 3.55 

3.3 0.9 0-50 

3.8 1.0 3.49 


2.6 3.3 630 

2.8 3.8 760 

2.8 3 9 800 

3.1 4.9 960 

3.5 6.1 1270 

3.7 7.4 1570 

4.6 10.6 2530 

5.2 13.7 3040 


Bigh Pressure Runs 

72 0.40 4.6 

76 0.35 4.5 

94 0.40 4.5 

77 0.28 4.5 

87 0.24 4.4 

103 0.22 4.3 

121 0.16 4.2 

110 0.12 4.4 


2.3 

1.0 

0.50 

2.5 

1.1 

- 

2.5 

1.0 

0.50 

2.7 

1.0 

0.43 

3.0 

0.9 

0.50 

3.3 

0.9 

0.50 

3.7 

1.0 

0.45 

4.1 

1.1 

0.50 
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8.0 CONCLUSIONS 

From a survey of .he li.ca.me, iiule work was found on high enthalpy supersonic 
parallel cen.ral injection flows, downstream of .he po.en.ial core. Of .he work 
reported, .nos, deah wUh i 6 ni,ion or reaction linuted flows. A detailed study was 
undertaken into ,he effec, of fuel injection parameters on ,he mixmg and con, bus, re in 
a cons, an, area due, with mixing limited high enlhajpy flows. Transverse ,nj 
experiments were also performed for a comparison. 

8.1 Displacement Th ickness 

The CFD parameter study showed that the wall pressure generated was mostly 
dependent on the stagnation temperature. This was because the temperature had a 
dlec, effec, on the fuel velocity which was the most important parameter ,n generating 
turbulence. The mixing rate was virtually independent of the stagnation pressure o e 
fuel as it has little effect on the fuel velocity. 

The displacement effec, of the mixing layer on the pressure rise was studied. It was 

found .ha, by reducing the area of the due, by the combined **»»>< 

of the mixing layer and boundiuy layers, tha, the pressure could be predre ed. The 

effec, of hea, addition was found to increase the displacement threkness by lowenng 

the density of the mixing layer. 

The effec, of mass addition through bounda* layer growth and mixing layer growth 
was found to be equally important to generating pressure* the hea, addrtron from 
ZZ ion. I. was found tha, an increnrenta. increase in displacement threkne* 
produced a different amount of pressure rise depending on how close the total 
displacement thickness was to the height of the duct. 


gj> Temp ^r 31 ' 1 ^ Effect 

From examination of the literature, i, was concluded tha, the effect of increasing the 
fuel stagnation temperature was its effec, on the fuel velocity for truxmg hnuted flow, 
and its effect on ignition for reaction limited flows. 

From dm results, only four runs had problems with ignition delay. These were all with 
cold injection. Heated injection a, the same free stream condtbons reduced the rgntbon 
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delay by a least an order of magnitude. This was consistent with the observations front 
Neer and Drewery (1975) and Huber et al [1979]. 

Increasing the fuel stagnation temperature increased the velocity of the fuel stream. 
This decreased the velocity difference between the fuel and air streams which in turn 
reduced the mixing. This was consistent with the wall injection experiments of Hyde et 
al [1990] 

The CFD results showed that the duct was very mixing controlled. The highest overall 
efficiency occurred with slow fuel injection, which increased the mixing efficiency and 
cold air conditions, which decreased the flame temperature and decreased the fraction 
of dissociated reaction products. 

The velocity and density profiles showed that as the fuel stagnation temperature 
increased, the velocity gradient decreased and the mixing rate decreased. The 
decreased mixing decreased the displacement thickness and resulted in a lower 
pressure rise. 

The different pressure rises between air conditions with the same amount of muting 
was explained by examining the displacement thickness. The different mass flow rates 
of air for the different conditions resulted in different areas under the pu profiles. For 
the same mixing layer width, this caused an increase in the displacement thickness with 
larger mass flow rates and a decrease in displacement thickness for lower flow rates. 

In conclusion, the highest duct exit pressure was with low enthalpy air and low 
enthalpy fuel. This condition was close to the ignition point. Increasing the fuel 
temperature decreased the ignition length but also decreased the mixing rate. 

8.3 Duct Height Effect 

It was concluded that if we know the relationships for the displacement thickness of a 
mixing or combusting jet and for the boundary layers in a duct we can calculate the 
pressure in a duct of any height with the same jet by simply reducing the area of the 
duct by the total displacement thickness. If the growth rate of the mixing layer and 
boundary layers is small so as to not produce strong shocks, then the pressure can be 
calculated by one dimensional isentropic theory. 
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It was concluded the combined effect of mixing layer displacement thickness and 
boundary layer displacement thickness that causes the change in pressure between 

ducts of different sizes. 


fi , 4 Equivalence R atio Effec l 

The pressure results show that the rate of wall pressure rise grows linearly to an 
equivalence ratio of 0.65 and then changes slope. The rate of pressure nsc then grows 
linearly from <>=0.65 to $=2.4. The rate of increase of displacement thickness with 
equivalence ratio of the air+fuel experimental and CFD results after $=0.65 was the 
same as the nitrogcn+fuel. This suggested that the further increase in wall pressure 
after $=0.65 was due to mass addition rather than heat addition. This was supported 
by the temperature of the duct exit which peaked at $=0.75. 

Increasing the mass flow rate of fuel decreased the amplitude of the experimentally 
observed pressure fluctuations. This was consistent with the observations of Gilreath 
and Schetz 11971] and Sullins and Anderson [1982]. 

Increasing the mass flow rate of fuel caused an increase in ignition delay. This was 
consistent with the observations of Synder et al [1965] and Huber et al [1979] 

The CFD mixing efficiency increased as mixture became more fuel rich or fuel lean. 
This was consistent with observations by Northam and Anderson! 1986]. 


§ c r^m parUnn of Trans yejsilgie ctipn and P ai a MInitcliQa 

Transverse injection experiments with cold and heated fuel injection showed little to no 
effect of fuel stagnation temperature on the duct exit pressure. The heated injection 
created a much greater disturbance to the oncoming flow resulting in large fluctuations 
in the measured wall pressure near the injector. This is consistent with the work of 
Chams and Collins [1970] who reported that the stagnation temperature had no effect 
on mixing but increased the radius of the bow shock formed by the jet. 

When the waves were normalised by dividing by the nitrogcn+fuel runs, the 
proportional pressure rise showed remarkable similarity between all fuel and free 
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suc am conditions. An ignition length of approximate,, ,50mm was recorded for aU 
conditions. 


When the transverse injection was compared to parcel injection with the same I fuel 
andTir mass flow rates, similar pressure rises were found. This was consistent with the 
simple mixing theory of Northam and Anderson [1986]. 


The longer ignition times of the transverse injection with cold walls 
consistent with the observations of Huber et al 11979]. 


and cold fuel were 
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6.1 Introduction 


For reasons outlined in Chapter I . supersonic combustion of hydrogen has received new 
interest in recent years. As a result, the amount of research literature devoted to the feasibility. 

„ .,*** .-*->» «,*. has ^ *. *«. ^ <~ 

books 'are dedicated to the review of current research. It is not feasible to go mto detatl on 
scramjc. design here and the reader is referred to Murthy and Curran (1991) for a review r of 
the topic The purpose of this chapter, then is to present concentration measurements taken 
with the mass spectrometer across the fuel jet downstream of a parallel injector ,n hypersonic 

Scramjets are essentially simple in geometry because the energy associated with 
hypersonic flow can bum through thin or intrusive elements of the engine. Because of the 
inherent simplicity of shaf*. scramjc, combos, ore can be modeledforerpenn^s,ud,esaa 
reduced scale and well approximated by rectangular ducts (Stalker and Morgan. 1984). 
Extracting optimum performance from the engine relies upon maintaining control over the 
mixing, and stabilixing the combustion in a short region upstream of the thnis, surfaces 
(Stalker. 1 99 1 ). The control of the mixing and combustion can pnncipally be achiev ed throug 
the means of fuel injection, making the study of jet mixing important for scramjel design. 


6.1.1 Parallel Injection from a Strut 

Four principle means of injecting fuel into the supersonic air flow have been proposed for 
scramjet engines: 

1 . Parallel injection from a strut or struts in the frce-stream; 

2. Parallel injection from a slot in the wall of the combustor chamber; 

3. •Vectored' injection from orifices flush mounted in the wall but inclined to the 
flow direction; and 

4 Transverse injection from orifices in the wall of the combustor. 

Of these methods, parallel injection from a strut is the simples, and the obvious starting 
point for a mass speclrometric analysis of supersonic combustion. Tins ,s because the model is 

two-dimensional and has the simplest geometry with no wall effects. 

Parallel injection from a strut involves passing hydrogen through the downstream end of 
a long narrow injector aligned parallel ,0 the air flow direction (Rgure 6.1). The two co- 
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(lowing warns .1.0 mix a, a ra,c .ha. is dependent on .he vcloo.es. densi.ies. compress, bill, y. 

initial turbulence of each gas stream, and the reaction rale of .he mixing gases. 

Feni Morel, i and Slu.sky (1965) firs, no.ed ,ha, if.be .cn.pcra.ores were hrgh enough 
, Uc „ as , he' fuel and ai, mix diffusively, .he ra.e of reac.ion was con, rolled by .he tale a, winch 
ibe diffusive mixing of ,he hydrogen and oxygen occurs. This ,s known as a d.ffus, 
controlled flame’. An al.crna.e form of.be diffusion com.ollcd flame in supe^mc flow occm 
when the fuel mixes a, lower tempera, ures and densi.ies so ,ha, the .gnmon ,s delayed. In .hi 
ease, after ignition. ,he diffusion of hea, through the premixed gases controls ,he rate of hea, 

release and the location of the reaction region. , . . 

The advantage of s.m, injection, as opposed ,0 parallel wall injection ,s that the fuel* 

introduced into the centre of the free-stream so tha, the fuel je, mixes on both stdes and he 
(lame front is kept away from the engine structure. Experiments have shown that where the 
(lame front intersects a wall there is much increased local hea, transfer, (Kanury. >• 
c Wots in shock tunnels have also shown that the combustion pressure nse from parallel 
ist er than parallel strut injection. (Stalker and Morgan. ,984,. The major 
disadvantage of the strut in the free-stream is the hea, transfer to the strut. Transverse 
injection from an orifice forms a complicated, three-dimensional flame structure which has a 
faster mixing rate bu, may cause too rapid a combustion a, low hypersomc n.ght speeds. 



Figure 6.1 Parallel Fuel Injection in a Supersonic Flow 


6.1.2 Theory of Supersonic Mixing and Combustion 

Understanding of the physics of supersonic turbulent mixing 
driven by empirical results. The equations of motion of compressible 


has so far been mostly 
turbulent shear layers are 
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milviical solution has been found. Numerical simulations are also 
50 complex .ha,, as ycl. no an > ^ ^ ^ jn Wghly comprcS sible 

currently unable ,o -.dcqua.e y r injection has been successful a, an 

srsrrsnL- — — ■ 

verity has led some marchers Orb* 

a c 063- Cassacio 1964; Zciberg and Kaplan. 1965) ,o attempt to aecoun, for the 

and Scheu. • transforming Ihe incompressible plane to ,he 

change ‘ ^ - ■“* *"* ^ 

compressible plane using hvDCr sonic flows even in the absence of chemical 

s T,r" 

reacnon . ( V- to lransport properties are. in themselves, important to the 

compressibility th $ 199)) , B rown and Roshko (1974) found that gas 

production of turbulent mix. g ' Ia ,e of the mixing region thickness 

co m p rc s s : ,i, r:,r^it:!oiams. « * — < - -*«* * 

srrru- — - “ fun T ex 7 men,ai rcscarch 

previously unexplored Row regimes is a valuable aid to further understandtng. 


6.2 Experiments and Results 

i nrder to examine the amount of mixing and reaction of hydrogen injected into a 
" • i am the mass spectrometer was used to measure the species concentrations 
hypersonic aus re . ^ of , hydrogen jet. For the same experimental 

a, several transverse local e was nKasu ,ed. ,0 compare the momentum 

Dux je, thmkness meler for these experiments was supenor ,0 any of the 

previous^ experiments, although the higher densities sometimes caused arcing between 
electrodes. 


6.2.1 The Experimental Test Flow Conditions 

Four different nominal test conditions were used. That is ,0 say that the initial shock 

, ub e^were^ = ™-— 

,n the Final cal “ la, ' d ^ order , 0 produc<! sufficient levels of static pressure for 

diaphragms in t s oc ^pies „f the four conditions are listed in Table 6.1 

— s in the test 

C T ^ ^5 TJZZ r f^tl: ■ “ I Mach number of 

:rrgh the combustor was Mach 55. the static tempera, urns of the flow ranged fmm 
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,000 ,0 HOOK, which was ,hc expected range of op.in.um .empennure in .he combos, ion 
cycle The densities chosen were ,he highes, obtainable a. which ,hc mass spectrometer couhi 
1 made ,0 operate, given the problems encountered with rite arcing and spread, ng of 

peaks These densities ranged around 0.05 kg/m . 

The expected contamination levels of ,hc test gas a. Ihcse condutons were, based o 
„„ of c ., aDlcr four, less than 10 * by mole fraction. The exception was the 

ITgZemhalpy condition which was expected to have a contamination level ° f “° und * 
by number. A, ,he highes, nominal enthalpy of 12.5 Ml/kg this may hate reduced .he AM 
sLiftc enthalpy to around 1 1.5 MJ/kg. based upon the room temperature spec, fie heats of 
helium argon, and air. As the influence of the driver gas could no. properly be mcorpora ed 
into the calculation of the nozzle now. the nominal conditions are used. As the mass of 
win is e 0 . small, the inHuence of the contanrination on the energy bounce m the now was 
I J,„e transport properiies could be altered. 1, was possible tha, the reduetton m 
“ ' Ires and pressums caused by the loss in enthalpy could be off-set by an tncreased 
:;r:: lema, Lrgy due to recombination in the nozzle ffow. was really only a concern 

a ,hc - « - - <*« p— * - — 

varied with time, so that in the interval of 0.5 to I msec after the star, of now, the stagnation 
pressure generally fell slightly. The amount of this variation was kept ,0 a mtnunum by he 
Ldition of argon to the driver gas (lacobs 1993). bu, the pressure was st.ll no. exacriy 
onstan. The variation can be seen in Ftgute 6.2 which shows the stagnatton and p„o, 



Figure 6.2 The Variation of ,h. Pltol Pr«sure a„d Stagnatlon Pressure 

with Time for a sample run at Condition D (Run If ) , . 

(The test time is defmed from 0.5 ,0 1 msec and the stagnation pressure has been delayed by 
' 0. 19 msec to allow for flow starting) 
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r ,;m e from the start of flow to the end of the defined test region of 
pressure as a funetton , o prC ssure ^ , hc pUot prc ssute is included in Table 6.1. 

now. The vanahon in the ^ , cvcl of thc now as ob,aincd from .he line of 

having been oblainc ^ Generally .he change in pressure level was less 

r The duration of .he .es, period corresponded approximately .o 3 passes ...rough ,he 
i - c ; n(J u pass was negligible in this instance. 
m0<,C, C NEWJcLlIns were ob.ained wi.hou. taking into account .he noale boundary 

, for the highest enthalpy condition, the possible con.amina.ton of .he test flow, 

layers or. for 6 P ? ^ ^ ^ ^ NENZ p calculalions probably over- 

According .0 .he result P ^ unitre! , im a,ed .he levels of nitric oxide. To 

es.ima.ed the £ o * chemisuy. .he norzle How calculations were also 
examine .he influence o conditions. These calculations 

performed with lhan were obtained in .he nonequilibrium 

gave a static tempera energy of the flow. Non-equtlibnum 

-> - - — * - — 

modeling in NENZF. 


Table 6J The Experiment al Test Flow Conditi ons . 


rarnr.mriTTi 
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The measured pito. pressure and srarie pressure for each of the eondWons are also go, 

in ji, r « r:r: -r^:ts:rrr 

measured pressurMwem' lower than rhe ealcu.a,ed pressures although .he fractional difference 
was within .he combined errors in Ihc wo values. 


6.2.2 The Experimental Model Configuration 


~ *» - ' - ■ ZZX'Z'Z r— T» ™ 2 

by ,64 mm high- Th, core had , racM ex.cn. which had no. properly been 

the duct was 91 mm. The deswners had a radial extentof 72 mm. 

quantified but based ^ (Jacob 199 „. The , oca, ion of the norrle 

83 mm downstream P ^ ^ ^ ^ 4Q ^ up5l(e am (,.e. 197 mm upstream of 

exit plane was , heItfore par. of .he .es. core of the norrle. The 

the injector exit). Some > hav e shown that the velocity 

- — - - - - - - - 

extremely because o. £ „ „ oid coroplicating probes with bounds* layers 
The large sire o c 199l) tv large duct also provided as 

which have occurred md^slu KS^ ^ Thc widlh represented die larges, injector 

close a simulation o J , unifoml flow, and the heigh, was chose, by 

which could reasonably be used P audy shock induced mixing 

Buttsworth (1994) to accommodate a shock pr S 

and ignition in high energy air flows. 
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As the ratio of the height of the duct to the duct width is 2: 1 , the duct cannot be said to 
be truly two-dimensional. To produce a two-dimensional free-jel in the shock tunnel test 
region was not a possibility however, and this configuration represented the best 
approximation. The jet may reasonably be considered as free in the transverse direction since 
the distance to the injector from the walls in this co-ordinate gave a line of influence traveling 
along a Mach line which intersected the centre of the duct at 370 mm downstream of the 
injector. This meant that the influence of the upper and lower walls would, at most, be a single 
weak compression wave from the boundary layer growth. Reducing the height of the duct 
would allow these waves to reflect across the duct several times in an equivalent distance. 
This was what, in fact, happened in experiments by Casey and he reported this as one of the 
major drawbacks of his experimental model. 

The strongest wave present in the duct would have been the wave originating from the 
front of the injector on the inclined side to the flow. This wave would have been able to reach 
the centre of the duct after reflecting off the upper wall in the distance beyond which the 
concentration measurements were taken. The effect of the injector therefore, was minimised 
by having only a short length for the growth of the boundary layer, whilst also avoiding the 
problems of the bow shock interfering with the mixing of the free-jet injection. The only 
influence which the injector strut provided was the slight mismatch in pressure in the two 
streams because of the different strengths of the shocks on either side. 

The injector strut, being 152 mm in length and 7.59 mm thick, is shown in Figure 6.4. 
For reasons to do with the experiments performed by Buttsworth, the injector strut was 
asymmetric, having an angle of 6.6 degrees on one side and 0 degrees on the other. The result 
of the asymmetry was to produce a slight mismatch in the static pressures on either side of the 
injector. This was calculated by Buttsworth as being of the order of 0.1% difference. The 
measurements were taken on the side of the injector which did not experience the bow shock 
off the 6.6 degree wedge. 



1 52 mm 


Figure 6.4 Diagram of the Hydrogen Injector 

Hydrogen was supplied to the injector from a Ludwieg tube 25 metres in length with an 
internal diameter of 25.4 mm. The hydrogen was released from the Ludwieg tube by a high 
speed valve actuated from the injection pressure. (Stalker and Morgan. 1988). The opening 
lime of the valve was about 1 msec. Measurements of the dynamic pressure between the valve 
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. that Dealt pressures were reached around 3 msec after t c v vc 
und .1* injector showed th « f* ^ ^ ^ onhc tcsl now and was triggered 

opened. The vaKe was o| ^ ^ ^ ^ m01i0 „ in the compression tithe. At the nine 

:“l.°n!w occurred, the pressures measured by the dynanne pressure transducer had 
fallen by 8% from the peak values^ d nozzle to a Mach number of about 3. 

Tlmrear cnds of the injector wer^^O. 19 nmrt ftom calibra iions performed by 

The conditions in the hy ^ ^ ^ cxi , p)ane of ,„e injector 

Buttsworth (1994). Bultswo of [cse „ oir rrfling pressures and obtained empirical 

a „a a, the throat of the injector for a range of res ^ ^ J ^ prc5 sure. This 

curves for the static pressur ^ ^ p[essure |osscs in , h e now from the Ludwieg 

empirical relation was ^ ^ „ ozzle was 4 ., 4. which determined the 

tube to the injector. r oogK gave a static temperature and 

Mach number and assuming a The theoretical Mach number from 

— “ — nts gave a Mach number of 3.0, 

». «. <* ** ***• ■<*** « ^ ,n ™ e 6 2 - 
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result of this nnsniakh was an exp of ^ effecl of lhe pressure oscillations 

of the hydrogen jet at the end o gcomctIy , after the initial perturbation 

propagated away frotn the ~ ^ d<wnsl « am . The need ,0 

Pl ayed no further inlluenee on h . — lhe conditions in the hydrogen jet 
adjust the state of tlte hy rogen ^ ^ ^ a „ isc „, ropic compression from the exit 

r sr : r- 

' u— - *■ - *• "" 

in Table 6.2. „ rPf < 500 mm downstream of the injector 

The profile of the injection wake was -asored 500^ ^ ^ 

which represents 66 injector heights do * n * U * a ^ ^ ^ 2Q ^ as measured from the centre 
of the concentration measurement were t ^ tha t the flow was not parallel 

of the duct. The distance from the ^ veloci ^ profl ie from the injector did not occur 

introduced the possibility that t e ce wa$ cxa mined by the pitot pressure 

exactly at the location of the centre of the duct. 

measurements. 


6.2.3 Measured Pitot Pressure Profiles 


thc velocity profiles across the mixing region with the concentration 
To compare the veto y P cocfflcicnt , the pitot pressure was measured 

profiles, and to obtain the value oft P achieved using a rake of in-line 

for each of the four experimental conditions. ^ cxlcma i diamcter and 2.04 

transducers connected behind a stainless steel tubing^ ^ ^ ^ ^ ^ ^ 

mm internal diameter and spac mm a pc B H2A piezoceramic type transducers. 

The transducers used for these measurement were calibralion was US ed 

which are standard instrumentation for s oc unn ^ ad j uslc d slighdy by a relative 

for the sensitivities of the individual gauges u ^ ^ app i icd t0 all gauges 

calibration technique (Stacey. 1989). The same pr measured pressure as a 

<— ' — — — 

function of time co-incided. T 8 . . . . . . fcw percent at most. 

the corrections to individual this mct hod as the error in the measurement 

Absolute calibrations were » ^ ^ ^ ^ of ^ absolule error in the 

of the calibration pressure was • causcd b thc expansion 

- *• — - - - 
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were in order of magnih.de Urge, than could be obuined in .be eabb,n,or. .he 

The signals prov.dcd by P P of , M|1; aIC sll0Wn j„ Figure 

Examples of .he onfiHered .he s.arl of no*. I. is likely lha. .he noise 

6.6. The ..me shown K r , ht n „ w> but instead .he resonan. oscilla.ion of 

*“ n °’ a -T “ Isducer The dimensions of .he .ransdoeer eavi.y were 2 nan 
lhe eaviry ahead of he Ira __ ^ 0 , which ,he inilial lenglh was 21 mm. Using 

diameter with a total long ^ ^ #f ^ 800 (he se dimensions gave resonan. 

,hccaleuU.ed fre«,ream^ ^ ^ ^ ^ or 2g 38 kH a. The observed dominant 

frequences of 300 <00 wh , ch COII esponded to .he lower predicted 

frequencies were approxim y This was also supported by the observation that the 

frequencies for longi.udina. '“-"Jh o.^c frn much of L ,es. period. Were .he 
pi,o. pressure ,n ^ n( genuine fie e-s.ream turbulence, .hen .he scale of 

Z" s,mel would have .0 be in .he order of a hundred millimetres ■ larger than .he 
dimensions of the mixing region. 



section A— A 




end view 


Figure 6.5 The Pitot Pressure Probe Rake 
(Frorn Butts worth, 1994, Fig 4.14) 
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Figure 6.6 Sample Time Histories of the Pitot Pressure in the Mixing Region Flow 
ngureo. f for Run 148, Condition A 

(Samples from locations 25 mm. 10 mm and 0 nun from the centre of the duct) 

The measured values of pitot pressure are shown in Figures 6.7 to 6.10 for the four 
different test conditions listed in Table 6. 1 . The measured pitot pressure at the entrance of 
due is shown as a reference line, as is the calculated value of the free-stream pr.o, pressure 
provided by the nonoquilibrium NENZF calculation. The predicted and measured values 
„i,e well bu, the pito, pressure a. the edge of the jet shows a consent nse to larger 
!,s 0f pitot pressure than the expected free-stream values. This observation is cons, stem 
1 previous measurements of pilot pressure in hypersonic ducts and md, cates that the free- 
stream gas has undergone a compression. This compression ma, have resulted from the waves 
1 (he boundary layers on the walls and from Use injector shut. Alternately, the compress, o 
in the free-stream rose from the entropy rise due to the mixing process and the subsequent 
gmwtht displacement thickness of the ndxing region. The magnitude of the tncrease was 

greatest general shapc of the profiles is that the measured values 

of pitot pressure at ±5 mm are not equal in three out of the four cases. The implication is that 
the axis of symmetry of the pitot pressure profile had been shifted slightly (- mm) into a 
positive location in the duct. This corresponds ,0 a shift towards the of the myemo 
Without the 6.6° wedge, which is away from the side of dm injector that has expenenced a 

stronger shock and therefore should have a slightly larger static pressure. 

The displayed error in the measurements represents the variation m the mean value of 
oitot pressure during the defined test period. As the pressure fell sbghdy dunng the test time. 
L mean value of the pressure was determined by finding the line of best fit through the data 
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using least squares regression. The standard deviation of the pitot pressure signals was much 
larger than this variation in the mean level. As discussed earlier, this large spread was not 
representative of the free-stream variation and therefore did not indicate that the error in the 
mean value was large, only that the oscillation of the pressure in the gauge cavity was quite 

large. 


Table 6.3 Measured Values of Pitot Pressure across Duct at 500 mm 


Downstream from the Injector 


“ — " 

Condition A 

Condition B 

. — 

Condition C 

Condition D 

Inlet 

601 

633 

623 

571 1 

NENZF 

625 

615 

605 

618 

Location (mm) 

. — 1 


-5 

351 ±56 

480 ± 44 

388 ± 29 

357 ± 14 

0 

307 ± 30 

431 ±21 

354 ± 13 

345 ± 04 

5 

323 ± 39 

433 ±78 

382 ± 22 

408 ±21 

•o 

408 ±61 

481 ±09 

412 ±33 

505 ± 13 

15 

607 ± 52 

602 ± 40 

519 ±46 

607 ±30 

20 

719 ±75 

654 ± 16 

617 ±69 

624 ±03 

25 

688 ± 84 

668 ± 33 

— 

— 



Transverse Distance (mm) 


Figure 6.7 Pitot Pressure Profile across the Mixing Region for Test Condition A. 

(Air Enthalpy 12,2MJ/kg) 
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Figure 6.8 


PHot Pressure Profile across the Mixing Region for Tesl Condition B. 

(Air Enthalpy 9.4 MJ/kg) 



X" . 0 9 Pitot Pressure Profile across the Mixing Region for Test Condition C. 
F,gure 6.9 Pitot Pressure ^ ^ ^ ? g MJ/kg) 
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Figure 


6.10 Pitot Pressure 


Profile across the Mixing Region 
(Air Enthalpy 5.6 MJ/kg) 


for Test Condition 


D. 


6.2.4 Measured Concentration Profiles 


.2.4 Mcasurru 

, , he measurement stations across the hydrogen jet were round by 
The concentrations at B£Cause only the concentrations during the test 

mass spectrometer rn Runs recolder was increased to allow more samples 

: were of inleies! the resolution o reducing ^ error enc0 unlered earlier with 

* <*» '" i,hin Umi " ^Tlw duration of the data recording was only 2.5 msec at the 
thin peaks such as the ium. . ^ Wng ftom 0 5 , 0 1 msec in any case. The 

her resolution bu, * d wilh mos , speclra. showing complete separation 

dity of the spectra was gc y Some spectra produced hydrogen peaks which 

the nitrogen, nitric oxide and oxygen !*• stages of the flow 

id the concentrations measured for these expe^ spcctromctt , measurements, a word of 

Before presenting the an ysr ^ ion g^nd ^ ejector were striedy 

aulion must be given. The con m wbcre ,he mass separation calibrations 

- “ 7 C °; 7 TTrJipal.y, the pitot pressure and the Mach number of the gas 

courted (see Section 3.7). Pnncp 1 free stream. Also, the wake region 

•xperirmced m^i^diarMte^l^h^frorrUhe pitot 

sss — - « - — 
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, iKBlicible with regard 10 <hc cff ' cl on roasS 

UPO or .6%. Hus ““ m now ins|dc , hc skimmer. (In regions closer ,0 .be foe. 

separation and the freezing properties can become very large in a transverse distance 

injector, where .he '" mK , sanl pling system would need serious attention to be 

of 2 nun. the use of the mas pc lhc size of the front skimmer), 

justifted. The problem could be r ^ ^ was only calibrated for the telattve 

Of the species of mtere . m ^ respect 10 nitrogen. The lack of 

concentrations of hydrogen, argon, ^ , ions had , 0 be made about the behaviour of 
comprehensive call bra, iondau mean could * nia<lc 0 f ,„e relative number 

— : - - - - - * -- 

molecules throughout the 'wn^sl^cTwareri the hydroxyl radical, and the free 

The remaining species qH, or H). The assumption was made 

atomic hydrogen, oxygen, and nitrogen ( , ■ • 1 ^ ^ ^ ^ calcu | a ,ed levels of 

that there was no free atomic spent* ze[0 in #, conditions. The calculated levels of 
free-slream atomic nitrogen using fo[ condition A. This is supported by the 

free-stream a-omte oxygen were s „ lhe siK «he atomic peaks relative to the 

results of Chapter 5 winch showed inconsistent vanat, on at 

molecular peaks as the enthalpy ^bustion reac ,ion rates predicted that 
higher enthalpies. Also, theore i ^ ^ , 9g2)> implying that the 

the reaction time was of the or er o ^ experimentally measured areas of the 

residence time of reactive species wa • ’ , we re not large compared with 

peaks produced b, ariiTteial depletion of these species 

the signal from diatomic speci . . was a largc contribution to the number 

relative to other species ^ was i® nQt discoun Ung the presence of such species, their 

rirs^rr^, fractions - ^ ^ 

beyond the error already present ."the* caUbtat ion factor for the measurement 

A u# assumption must ^ as ^hing else, this was taken to be 

of the water molecule. For reas P* ^ ^ B , he calibration constant for 

The actual calibration cons,® ^ itot pressures and the calibration constant 

helium (the larges, measured) a.) ^ ^ monatomic and the mass ratios for these 

for hydrogen is I also. toW* ^ ^ nilroe e„. the error in the esrimate of I 

species are more severe that ^ 25% w fcch j s to say that the calibrauon 

for the relative c#ib,a,ion of assump,ion the nmle free, ions of 

might be as much as E25 or as ^ wividoal rclativ e concentrations to niuogen. 

the seven species can c cu calculated from the areas under the peaks in the 

The fractions for the seven species were caicu.a 

recorded spectra and used the following formulas: 
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(HjI 

IN,) 


1,68 Area jji 

1 x 0.557 Area N, 


(6.1) 


to 2 l 

IN,] 


1.68 Area °2 
' x 1.57 Area N, 


( 6 . 2 ) 


[NO] 

IN,] 


1.68 Area NO 
1 x 1.89 AreaN, 


(6.3) 


(H,01 

[N,l 


1.68 Area H 2° 
1 x 1.31 Area N, 


(6.4) 


Driver 

IN,) 


1.68 Area He 
17 x 0.313 AreaN, 


1.68 Area Ar 
+ 1.25 x j g -j ^eaN, 


(6.5) 


. , ,, mo , e fraclions are obtained by dividing the relative mole fractions by 

Here the absolute mole fracltons ^ ^ f()1 , ^ ot , he 

; sum of the relative mole fractions o a 1 1 to 6.13. In these figures the traces are 

ns as a function of time t are ^presented m of ^ (b) and the driver gas 

vided into hydrogen and P traces show a variation 

MmS ec after the stan oftheflo^ region „ ave ,**„ obtained from the 

Concentration P‘°« ^ ^ „ show „ m Figure 6.14 to 6.17 for the test 

neasurements from £ , staniald errors of the mean of the measurements 

ondittons A to D. The err inc|ude systematic errors from the assumptions 

rom all spectra for each m f res is included in Table 6.5 The 

egarding the calibrations. The . R 6 .18 t0 6.21. In Figure 6.22 the 
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Table 6 4 The Condi.ions for .he Tea. Runs fur .he Measurement of Species 
iaDe ° ronccntration — 



B 


k±69±0.05K) 


10 


Uss + finffi I2t0.0)k)22t0.01 


15 


20 


10 

15 


20 


|0.50 1 0.03|0 


loos ± o.osto 


|0.68 ±0.18 


0.18 ±0.02(0 


lO66±0.04 


0 51 ±0 05 
0.35 ±^05 


low ±0.050 
030 


0 05 ± 0 
Io05~±0. 


k).13±0.l 


In 64 t 0.0sk).04 ± Q.ook) 23 ± 0j> 


0.47 ± 0.05 0.03 ± 0 00 0 35 ±0.01 


U ia ± o OSiO 02 ± 0.0d0.44 ± 0.01 


15 

20 


004 * 0-00 0 .47 ± 0.02 0.04 ±00 1 0^1^001 


10 


20 


.04 0.04 ± 0.06 0.46 ±001 0.06 ±0.00 0 15 ± 0 01 


0,68 ±004 
0.40 ±0. 


06 ± 0.01 0.18 ±0.01 0 03 ±0001003 ±0-00 


In 03 + 0 Qolo.02 ± O.QQlO.05 ± 0.01 1 J 


12 ± 0.0210.28 ±001 


07 ±0.0110.48 ±001 


.020 


0.04 ±0.01 


ln02 ± 0.04|0-Q4 ± 002 


0.01 ±0.00 


17 ±0.030 


22 ± O.OllO 02 ± 0.00| 


30 ±0.01 


0.06 ±000 


34 ±0 01 


n 06 ± 0.03 0.05 ± 0 00 0 60 ± 0 03 003±0jH 


U.HU i. U.uj — — ~ _ _ . 

0.24 ±cT 05 0.05 ± 0 02 0.49 ±0.01 0 0S ±001 OMdbOOi 


|0.70 ±0.01 


5.04 ±0.01 


i.06 ± O.OllO. 12 ± 0.001 


03 ± 0^0 110.05 ± 0 02 


1005 ±0,00] 


[0.08 ± 0.00| 


0.0910.01 


Uq U.W J.W WI v- . * _ 1 — ~ , Al 

ns n 05 ± 0.02 0.39 ± 0.0 1 Q OS ± 0.00|0.1S± 0 0 0 0 0° 1 0 °l 


!0.04 ± 0.02 


0.20 ±0.01 


0 05 ±0.01 


0.01 ±0.00 


0 03 ±0.01 


0.16+0 01 


lk02± 0001006 ± 000| 


lnos±ooolo.io±oool 


too7±o.ooloT3±ooo| 


K±16±0.02 


ft ->1 + n nn n 03 ± o 0d006 ± 0(XM0l ±000 


b.04 ±0.02)0.19 ±0.01 


0 04 ±0 01 


0.14 ±0.01 


0.01 ± 0.00 


000 ±0.00 


0.01 ±0.01 


^06±0.01 


0.11 ±0.02 


iooo±ooo 


io.03 ±0.01 


JO 

10 


to 


10 


10 


10 


10 


0.52 


1.00 


0.86 


1 21 


1.02 

1.48 


1.17 


1.76 


0.71 

0.51 


0.97 

099 


0.50 


098 


1.03 


1.31 


053^ 

0.51 


0.53 


1.01 

1.00 

1.01 


0.92 


1.08 


1.40 


1.56 


0.98 


0.49 


1.23 

0.97 


0.49 


1.51 


0.97 

1.73 
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Figure 6.11 (a) The Mole Fraction of Hydrogen and Water as a Function of Time 

(Condition D and 5 mm off Axis) 
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Fieure 6.11 (b) The Mole Fraction of Nitrogen Oxygen and Nitric Oxide as a 

Function of Time 
(Condition D and 5 mm off axis) 
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Figure 6.1 1 (c) The Contamination Level of the Test Gas 
(Condition D and 5 mm off axis) 
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Figure 6.12 (a, The Mole Fraction of Hydrogen and Water as . Function of Time 
* g (Condition A and 10 mm off axis) 
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Figure 6.12 (c) The Contamination Level of the Test Gas 
1 8 (Condition A and 10 mm off axis) 


72 



Mote Fraction 


SCRAMJET COMBUSTION SURVEYS 257 


Run no 120 


0.9 - * 

0.8 : i 


IS 1 


Hydrogen 


r r r : i 


).i r 

n 


Water 

D □ d a fl Oo nOOn n oofl a aa °o 0 °.° Q £Lo n-fl-riJ3- 


Timc after Start of Flow (msec) 


Figure 6 13 (a) The Mole Fraction of Hydrogen and Water as a Function of Time 
8 ’ (Condition B and 1 5 mm off axis) 
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Figure 6.13 (b) The Mole Fraction of Nitrogen Oxygen and Nitric Oxide 

as a Function of Time 
(Condition B and 1 5 mm off axis) 
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Figure 6.13 (c) The Mole Fractions of the Driver Gas 
(Condition and 15 mm off axis) 
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p[.ure 6.14 Mole fractions as a Function of Transverse Location 
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Fieure 6.15 Mole Fractions as a Function of Transverse Location 
8 Condition B (9.4MJ/kg) 



pi „ U rc 6.16 Mole Fractions as a Function of Transverse Location 
8 Condition C (7.8 MJ/kg) 


75 



Mole Fraction 


60 CHAPTER 6 


0 -* T Jh dr open 



Water 


# j|L=*a. 

» r a A 


Oxygen 



Nitrogen 



Nitric 0%ide 




Driwr Gas 


\ 




5 10 1520 0 5 10 


° 0 5 10 15 20 0 5 10 15 20 0 


1*520 0 5 10 15 20 0 5 10 15 20 


Transverse Location (mm) 

Function of Transverse Location 



76 



Ma« Fraction 


SCRAMJET COMBUSTION SURVEYS 261 



Trans wrsc Location (mm) 

Figure 6.19 Mass Fractions Condition B (9.4 MJ/kg) 



Figure 6.20 Mass Fractions Condition C (7.8 MJ/kg) 
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Figure 6.21 Mass Fractions Condition D (5.6 MJ/kg) 



Figure 6.22 Average Molecular Weight across the Fuel Jet Mixing Region 
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6.3 Discussion 

Generally the accuracy of the measurements was not as good as those of the driver gas. 
This was due to the problems with calibrations and performance of the instrument. The 
systematic error in the measurement of hydrogen was 20%. water 25%. nitrogen ,0% nitric 
oxide 10%. oxygen 10%, helium 5% and argon 10%. These estimates were dommared by he 
estimated error in the calibration constant. The random error in each measurement (due to the 
sampling rate or lack of peak definition) is included as an error bar for each measurement. 
Considering that this was the first lime that measurements of thts type had been performed , is 
kvei of error was no. unexpected. Despite this, there is still a considerable amount of 

information contained in the measurements, much of il unique. 

An additional difficulty with ,he results can be seen in the traces of the concenlralions IK 
a function of time (Figure 6.11. 6.12 and 6.13). That is. tha, the concentrations a, eac 
location in the flow were no. steady with time. The general trend was for amcomponeo. 
concentrations to fall with time after peaking a, around 0.5 msec of (low while the hydrogen 
conceniration rose with time after 0.5 msec. Where there was noticeable water the waer 
concentration fell with increasing time. An obvious explanation for this was that the effective 
heigh, of the hydrogen fuel jet (after adjusting to the ambient static pressure) 
with time as the static pressure of the test gas fell. Such variations were ."dependent of he 
operation of the mass spectrometer, and are therefore no. so much a P^lern w. t h *e 
measurement, as with the interpretation of the data. In fact, this observation htghl, glued the 
importance of the instrument’s ability to obtain time resolved measurements, as single instant 
concentration measurements would no. reveal the unsteadiness of the combustion. 

A comparison of the results from the different test conditions, tndicated that the highest 
enthalpy air (low. Condition A. showed a considerable amount of water across the fuel jet. A 
IL, Lun, of water appeared ,0 have been produced in Condition B.whtle m dte two owes 
enthalpy conditions, there was no measurable amount of water produced. In all thee of the 
lower enthalpy conditions there was significant penetration of the oxygen into the fuel je 
without reaction. In Condition A. the amount of oxygen measured inside the fuel jet was 
significant This indicated that a. only the highest enthalpy condition had complete reaction of 
the mixed air and fuel occurred, and that the combustion of the fuel was mixing-limned. In the 

.. . / R r and D) the fuel and air mixed but did not react, 
three lower enthalpy conditions (B, C. and U), n. 

Condition B may be reaction-limited, since there was some water present, whtle Condemn C 
and D showed no evidence that ignition occurred. In all conditions, the measuremcn of tog 
amounts of hydrogen within the je. indicated that nothing close ,0 complete mnung of the fuel 

,e. occurred. |k( ^ |he „ ansvcrs e extent of the mixing region taken from the 

concentration measurements and the thickness of the momentum flux defect both increased 
with the increase in test flow enthalpy. Although the initial effective fuel pit thickness also 
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increased with increasing enthalpy. ^dZZk fraction of the 

jc , height 5 “ ^ a ^ 1 " o ; ;“ iimlta ^s, y for .he four test conditions in Figure 613. To 
components of tur are plo frec . slream molecules into the fuel stream, the water 

obtain a measure of the pent. tj and oxygcn . Therefore the mole 

mole fraction of driver gas. 



100 \ 


n L— — - r- 1 — — 

•5 0 5 10 IS 2025 -5 0 5 10.55025 -5 0 5 .01520 


t-- ^ 0 

•S 0 5 101520 


1U last'*-' - 

- nktanee from .he Centre of .be Due. (mm) 


. w een the Pilot Pressure and the Air Mole Fraction 500 mm 

Hgure 6.23 Comparison b«l.«^eW» oriheInje£(or 

(□ - Pi.o. pressure; n - Mole fraction of tes. gas (including water and driver contamination)) 


r , am D i,ot pressure and the hydrogen jet pi.o. pressure were predicted .0 
The free-sueam pi ^ iijfkcKd m , h e pilot pressure measurements. 

increase as the A Condition A. the centre line pi.o. pressure 

except in the centre oft ) and rep[CS( . nttd , grca rer deficit in momentum (lux 

was significantly below It ex(* ^ ^ nux de( - |cU of lhe fad jet. 

(below the free-stream mom'" 1 ^ irldicatcd ,ha. the density was falling 
That this occurred coonc, ^ increastd U .fas was so. then the thickness of 

on the centre of the fuel jet a F . and the free-stream pitot 

- « * generated b, the expand of .he 

r**- « ■ «* ° f —• " c ° mpressi °” waves 
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- — flux profiles - ,e — £££ ^ 

roUSh,y ' qUiValC h "' rr,ull“ ion doming ovc, .he molecul* transport. In a free 
mixing region where lh= Kmi , nfinitt n „ids. .he grow,h ra,e of .he heigh, of .he 

turbulent shear layer ^ |incar , y ^ downslream distance. (Figure 6.24). 

mixing region has be ' £ ing dcn0 ,ed .he spreading coefficient. (Schell. 1980). 

with .he constant of pr P ^ now has ^ found l0 be a f„„«ion of the velocity 

and air streams, the density ratio, ,1m level of density fluctuations, the 
rand between ,1m hre in lhe shea, layer, the rate of hea, release in a 

b ri rl the initia. turbulence levels in the free-stream flows. (Dimotaltis. 
chenacally reacting slrcam ve|ocitics , where intuition would suggest that the 

ndxhtg would teshght. the initial turbulence and densiry Ouc, nations ac, to produce stronger 
mixing than molecular diffusion (Brown and Roshko. 1974). 


Air Stream 



s 

Figure 6.24 Definition of the Spreading Coefficient “ 


- -r. - “-err 

."-.TS— ™ - • •- w » ““ “ “ "" 

correlated to the function (Dimoiakis, 1991). 


^r,s) 


C ‘ 2(1 +s r) 


1 - 


11 


s , fl )/(l+s m ) 


1 + 2.9(1 + r)/( 1 - r) 


;) 


(6.6) 


. • .• c k the density ratio between the two streams. The constant 

ol";l'd unknown factors and has fallen in the range 0.2 5 <c<0.4, 
vTues of density and velocity ratios are listed inTahle 6.6 for test conditions A to D. 
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_ . , , , Thp Enmlrical Predictions of Spreading Coefficient and Mixing 
Table 6.6 The tmpirica. K ,„ rimc „,al Conditions , 

egion Height 


r.jc — 

Condition A 

Condition B 

Condition C 

Condition D 

Velocity Ratio, Uj _ 

0.5 

1.1 

0.6 

0.8 

065 

054 

0.74 

0.35 

Density Kano, Pj p» 

Incompressible SpreadinjCocfnc^m 

Oil 

1.2 

009 

1 

0.07 

0.87 

005 

0.6 

Convective Macn Nurnpo 

Predicted Compressible Spreading 

Coefficient — 

Pra.a-.ntd Mixing Region Height (mm)_ 

0.02 

11 

0.02 

10 

0.02 

10 

0.02 

12 


, n nn rpvtine shear layers has been suggested by 
^ T cy n reduced he concept of a convive Mach number 

irrrr if i no. - . - -» * - — - 

the convective Mach number is defined by 


U, + U 2 


Me 0| + fl 2 


(6.7) 


. „ M i* ihe conveclive Mach number, U is velocity and a the speed of sound, with land* 
where M t .a the co __ scn ts Mach numb er „f the large scale turbulenr 

denoung each * • For ,he experimental condUiohS of .he current 

structure with r J S numb e r i; iven m Table 6.6. Using bus as the defining 

measurements, the convect, ^ ^ , ha| by action 6.7 

rrrnvLTve K Mach number increases. The experimental measurements are correlated by 

an empirical relation: 


AM,,) = 0.8e' 3M " +0.2 


( 6 . 8 ) 


Where / represents the reduction in the compressible spreading rate from the 
...f S rc , ation is plotted in Figure 6.25 along with the experimental data which 
incompressible. This re P As can be seen the majority of data is taken at 

i, is mean, to represent (sources rediced from dre 

convective Mach num*rs lets £ are also given .Table 6.6. AscanM 
convecuve mach num Nations and the convective Mach number 

ZZ7SZ iivaJ predicted spreading coeffiCn, for each experimental 

condition. 
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• J hr re immediately behind the injector, two shear 

— :r:x r- *. ^ * ** « ,o ^ ** 

" T," « a .hou S h not necessarily independently from .be influence of .be other 
downstream distance. S WJWS |hcy produced. Each sel of waves would 

The two shear layers wou ““ iHy generating additional turbulence and 

travel through the ^ pom, downstream, the two shear layers would intersect 

increasing the ^ ctoge from dual she ar layers to a jedwake flow. After this 

and the nature bcfofC duc t0 lhc coupling), the empirical results from semi- 

had occurred. (thoug pe p • w (he growlh of lhc mixing region, 

infinite shear layer of , * , w0 sbear Uy ers. the properties on the centre of the 

Upstream o properties and the pitot pressure and concentration 

* """ "IT" t “ial values equivalent to those a. the start of the jet. This region of near 
profiles should show ^ ^ , corc . Dow „s,ream of the potential core, the 

constant properties ts r fre e. slr eam properties, and the decrease in difference in 

axial flow properties app distance. In the far downstream, the 

quantities usually '“''^^^m^ralure become small and the growth of the shear 

transverse r**"» ^compressible growth rates (Bradshaw. 1977). 

could be approximated quite wet y significant, the length 

- - * - - ■— 

of the potential core experiments, this length is 120 to 130 mm. 

'ZZSX2 ,hc i " icc ' or - 1 ,he "" n,s - ,he 



(Dimotakis. 1991. Fig 4) 
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mixing region will be behaving as a jel/wake flow and the empirical correlations for semi- 
infinite shear layers will not apply. It remains of interest, however to discover how much 
greater the growth rate is in a jel/wake mixing layer than in a free shear layer, and if the 
predicted independence of the spreading rate with experimental conditions also applies to the 
jet/wake flow. 

From the measurements, the half-height of the total mixing region was approximately 20 
to 25 mm for Conditions A, B, and C and approximately 20 mm in Condition D. If the initial 
distance required for the isentropic compression of the fuel jet to the ambient pressure was 
ignored, then the spreading rate was about 40/500 or 0.08 (0.07 for Condition D). This was 
much larger than the predicted spreading coefficient for free shear layers and represented a 
much faster mixing and entrainment of the jet and frce-slream gas than did a single shear layer. 
The incompressible spreading coefficient for free shear layers was larger than, but much closer 
to, the measured spreading rates. Tlus was consistent with the conclusions reached by Casey 
(1991) from strut injection measurements in a narrow duct. His experiments were also 
conducted far enough downstream for the flow to become jel/wake dominated. The 
independence with test conditions does not seem to apply to the jet/wake flow since the height 
of the mixing layer is seen to grow with increased enthalpy. 

Few other measurements of species profiles in shear layers have been reported, and none 
in a high temperature, high Mach number flow. The few measurements that have been taken 
imply that the transport of momentum, energy and species between the two streams docs not 
take place at equal rates. Chriss (1968) found that the centre line variations of hydrogen 
concentration and momentum along the axis of a free-jet followed different functional forms. 
Gincvskii (1966) found that the transport of heat occurred twice as fast as the transport of 
momentum across shear layers in incompressible flows. The current results indicate that for 
these experimental conditions the transport properties are very similar. 

Past research has shown that turbulence levels in the facility have an influence on the 
growth rate of the mixing layer. This was not quantified in these experiments, although the 
noise level in the pressure measurements was quite large. 

In summary then, the state of knowledge of the mixing between supersonic fuel and air 
behind a parallel injector is not sufficient to theoretically explain the existing measurements nor 
to predict the outcomes of new experiments. In such a situation, measurements of 
concentration profiles across the mixing region downstream of such a combustor must be a 
valuable addition to the measurements of density and pressure which are routinely obtained. 
Confidence in numerical simulations relies upon the ability of the equations used to model the 
flow being able to reproduce experimentally obtained results. For these reasons also, the 
exploration of new measurements in new flow conditions is essential to the development of 
better understanding and modeling of the physically dominant processes in supersonic 
combustion. 

These experiments in hypersonic combustion have only scratched the surface of the 
potential mass spectrometry offers to this field. The presentation of these results and their 
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■iv hrirf as further detail is beyond the scope of this thesis, 
consequences J ^ cxpeli „,ents pose more questions than answers and 

The information obtained from these ^ research is required 10 more clearly 

highlight many further » vc " u ” f °' '“'“lading rate in parallel injection in a supersonic flow 
define the mechamsms control', g JP . jn ^ , ffl delermi „ c whether the spreading 

as well as to sample the gas c „ w „uld also be valuable to perform 

~ ,o edn a 8rea,tr 

understanding of the turbulent mixing present in supersomc combushon. 
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ABSTRACT 

IT ,h, ,«»Us of an experimental investigation 

This p*p*r presenls ^ ' tuition to turbulence 

Into the characteristic* of bo **> experiments was conducted 
in hypcrvetooty uiD _,j with static pressure and thin 

using a flat plate mode ,tq PP” shock tunnel. Tnn- 

61m beat transfer ™ fin , e 0 f 2 35 to 

sition was “? Re^oW* number eorrelate* weU with 

192 MJ/kg The tran..t^Reyr»M‘ relation. 

— **** cooting are 

Errs srt 


NOTATION 


, speed of sound 

wave amplitude in unstable tegmo 
i wave amplitude »t tbe neutral point ,. 

[* amplitude of tbc sound radiated from the nozale wall 

1 boundary layer 

n constant of proportionality 

p diameter of tbe nozzle 

£■ normalised po*er density 

f frequency (Hr) 

/ peak frequency (Hz) 

t dimensionless frequency -fjv 

> dimensionless frequency of the most amphfied mode 

if specific enthalpy (J^g) 

L model length (m) 

M Mach number 

p pressure (Pa) 

p* atmospheric pressure 

p pilot pressure (Pa) 

q beat flu* (W/m J ) 

M^^ved HSS 30 ^ 
bee 1993, accepted 23 December 1993. 

Paper No 1939. 


<> 


specific gas constant (kJ/kg K) . . 

Reynolds number based on the distance from tbe lead- 
ing edge 

transition end Reynolds number 
transition onset Reynolds number 
unit Reynolds number (Wm) 

A M 

Stanton number * pt/(Ji, - O 
time (s) 

temperature (K) 
flow velocity (m/s) 

ratio of sped6c beats; interautieocy factor 
fraction of transition zone de6ned by Equation (9) 
driver gas compression ratio 
density (kg/m*) 
test time (s) 

kinematic viscosity (m/s) 
average quantity 


Subscripts 

0 stagnation quantity 

quantity at tbe wall 
r recovery quantity 

Superscript 

fluctuating quantity 


INTRODUCTION 


The prediction of boundary layer transition u l«««Uy . £*** 
related to high Reynolds number flows. Transition bar been 
ly investigated for many flow situatK.nt tn the past, but there has 
been link investigation into tbe transition phenomena relating to 
hvncrsonic high enthalpy flow*. The boundary layer has a l.gn.^ 
2 ? influenced the c”mal flow field, of re-entry brfj^ueh 
as the Space Sbuttk. and would have i critical influence on 
bypersonicduct flows such as ttose created in scramjet engines. 
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0 * cere of ,ub tonic and I ! 

boundary layer instability and u oumb «T and tbe I b 

S io the unit Rc^W, ^ o M»^^ wivel ^ 
cooling of tbe wall™. The * u, lt^ bv ?S . U bili«y theory « 

and revealed by “P*"®***^ f^Tunear bounds lay« 
transilioo i, narked by (be ^^^J^anc* of local < 
characteristic* which were “^^-^^Reynold, num- 
Ubuleoi (noo-Unear) eve "“_ 1 ^^^ w oveI flat plate* and I 

* '“rL^^ V^rrclated with the influ-U 
cooes, was of the orocr i\r. . 1VI .bowed an I 

parameier*. bw*i^oo jJ^^g Reynold* number. Mach 

increasing trend with uicre thought at one tune I 

Dumber And boundary V in complete boundary I 

i - ■ssss: ^ &» ■>' — *• 

I , 1^w-» W( g^JfSKS£t 

mental change* occur u ‘ faltaWBv wave*™. Lin- 

1 associated with the ! Je '’ e,op . addition to some low I 

1 ear liability analysis dcm 0 r f. ... , ves w bicb are dominant in 1 
frequency. or firs, mode ^re «sU« higher 

I subsonic and supersonic -stabilities, in hypersonic boundary I 
I frequency, or second mod • 6 rs i mode and I 

11>w ^v"- 

I becomes dominant in the process * cooling, as deter- I 

baity theory * 1 indicate, that enhance sec 

I mined from waU to stagnation Uyer and thus I 

I ond mode instability »n by(W«»f stability exper- I 

I reduce the effect of 6 rst mode discovered that the 

iments conducted by S ‘ e “^^^e,t,, equal to twice 
I wavelength of the seoon assumed that the wave velocity 1 

I the boundary ^ 0 * boundary layer edge Stetson 

I ^ KirruneT speculated that 

I lately 10 tado^Mtja lbeJ^J.^ j^hT hypersonic boundary 
I face cooling should delay distuibances were known to 

I Uyer. However, when sew ^^mode transition of 

I be dominant, cooling aupported by Seddougui’s™ ou- j 

I the boundary Uyer This ^ ^^«»lin* enhances the 

merical study which suted ^ accelerate* the 

Tollmein-Schlicbting viscous ^^XrUyer. 

MM iUonpmc«, in Ok compressible b^^Uy fe . 

Previous hypereon* transition ^nmenu ^ ^ Muf . 

I searchers have yseW ** 1 w4i little influence of 

ptyn. mid Ross* have .nd.«^t^ of ^ 

I waO to recovery tempera^ ratwJV,^ experiments conduct- 
I transition point. ^ling definitely caused a sub- 

I ed with a conical model that wan V. ^ Shock tunnel 
I stantial upstream movement of dbe tr demonstrated a reduction 

'V***™ H s«w« «r T/r. -» 

I is transition Reynolds num SioUaV H) showed that ibe 

reduced. The results of , nd re-revered, 

increase of TJT, may 4£si«e » ■» ° f 

number with wall cooling andne > tiper iment* 

All of the above ?nd hence 

I bad relatively low sugna , ( -p* , M , conditions of 

| tbe flow velocities were also relatively tow. i « 


t2 

conditions. 


equipment and ex perimental conditions 

Esoeriments were conducted at tbe Unirersity of Queensland 
ushT«X T4 free piston shock tunnel (Fig. I). A contoured noni 
s\ fi,,d area ratio was used and the nominal Mach number 

r^SOAn“pl'“o*“ (Fig )) «» . W *«» -*>■••• 

plac L in the test section. The blunlness of tbe 

was 008 mm. tbe toul length of the plate was 

width 230 mm. A fence (not shown in Fig. 2) was attached to e*4 

tide of the model to curtail the disturbance* generated at the 

^•fLTSU rff n- «»> 

r,om Ik. kUing «lg«. *c «t°«l »«•«*»• ■“ lM “ 

Consecutive transducers were 1 7 5 mm apart. 

C.u,-.k.« ..F..IPF oMke •» « *' '7'”S 

MJAg The nominal conditions were T» 

kPa and M = 5 2 - 6 *. The flow speed vaned from 21 

'.SUM* i" «* (>“•"" - **st2 'S 

f,« Ser»tures during the test*. *od temperature trace* 
then inte crated to obtain beat transfer rilea A pressure transducer 
wJ SSSled in tbe plate to obtain a reference static press-e 
and to determine the steady flow condition! _ 

-ft* main objective of the present mvesUg.Uon wu m de« 
th* transition phenomenon in bypersomc and high enlb ** 
SiSSrAowi. However, tbe enthalpy range over whwt 
« observed is limited by the lengths of detest core- 
thVmodel In shock tunnel tesU with a fixed nozrie. bigb eolhUpy 

is'achieved by increasing sugn-ion , «4 ^ **%££* 

KsgasaxL » *‘°';“** 

inoreased. the unit Reynolds cumber decreased and tbe Reynott 
^^Tbased on the distance from the leading edge dm 




Fi 0 ur# 1 . Sketch ol the T4 shock tunnal. 


Figure 2. Sketch of the Hat plate model. 
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decreased. TW«. ! 

downstream from the k«(W mgn8 tjon enthalpy of I ( 

beyond the model kn « MJ P* bul *, onset of iransilioo could be 
the How tested wu 26 MJK* «>ui w I 

deleted for stagnation en^slpw oy ^ # when a 1 

Tbe decay of tU 8 nlU( *\>*“. t ui ,«| for the examination of 

steady flow ov« * «"“£ 33-2 Is U u* ** 

U?; tts «*— * "° w to ^ tibe enure 

model U I 

i=£ ...(»> 

I Tbe behaviour of ***"*** ^^^ 0 "^^ m«Jel For the I 

“ten the flow .peed U the lowest, we have 

I ,-±1 = 2 9xl0"*s = »0l‘* 

I Over this period the rtlen ^^^* n fl oe ^ h^wnsfer mea- | 
I amount of decay did not su aj deeded as discussed later. I 

I turement through which reeion of the boundary layer I 

I It should be u^th*Tbe pressure decay during the 

•ZZJfcZi rsf r-*”* 

I To investigate tbe influence 20% P „ fcC n and 80% 

1 sition. a series of runs the fiow cond,. ions at I 

helium driver g« ««£ ™^*cay in stagnation pres- 
I an enthalpy of 3 5 MJfcg. ™ ^ detected Howev- 

I sure. No difference in transi , ^ contamination free 

er. bod. die maximum stagn-mo 
I test time achievable by mixed dnver gases wer 

I those achieved by helium *»»»■ , hock tunnel was 

I Tbe test lime in a hilbenlhalpy renec „ 

I limited by the driver gas ^ Ufcated reflected 

the vicinity of d* .bock «ubew»n undergoe** ^ ^ (he 

| shock, and suffers a smaller clungem sbock „ 

I centralised gas which P*“" * ^ te$l region along the 

1 therefore forms a jet ThU problem bad been 

zsrsss 1 ssT^TJ*-’- 


calculated the lest lime in lean, of test slug length 

C*JCUJW~ »« diiphrttm rupture pressise of 57 MPt 

.ugnauon er.^ Mlftg d* •“» »>“« w “ ,boo< 

<B * \irh was more than three timer the model length. Tbe »U| 

a ri sp- 

£££?* «-*>«* ^ • m %£,* 2 ?szLi 

r-s. *«* 

SSWl »■»«««■ *“ i “* 

Tbe test »cc , lf) Given the stagnation conditions in 

»*. a. n- 

the noizle r , , . as through tbe nozzle 

rK.?« 3E£Si.-S *— oic-g^ 

™ l»Zun>S for die flow upstream of tbe duoat.^mKd 
^ TK^ttm calculations began downstream of tbe throat and 

U^'cbemical kinetics provided by * n \^ J V ^* e ‘^ U ed 

^ for 'temperatures greater than 5000K. 

?z££?£zsr. K %=:r : s 

■ n J S5» j£e enthalpy rang* within which tbe transition expen- 
ment was conducted. 


FREESTREAM DIST URBANCE MEASUREMENT 

-»* "?jZ°s£L £J 4 S<%*- A KB 

i t^crZisES «. £.»* .*»* — rg- 

m ol 500 IH. T>* ■"»/»” 
ibe recorded data was not less than 250 kHt. of |£1| 

T £ variations or the fluctuation 

Pissfa 

TCS - -sr-SS.'i wts 

Lr'.S'i.xrrs «-r 

lion intensity of the velocity and that of pitot pressure is discuss*, 
below. 


AW**** iJourn * 


ra* 






ai ' 

V* unit R^°«* n^br. 

*- ** SrSfSSS. 

the pilot pressure r , u> * n >P eriu ‘ 
assuming the gas “ P^^ 1, 

/'IV 


P, = Cpu 5 


where C is a constant and 




ri Li 


Therefore, we have 


5l = (M + 2)^ 



f (kHz) 

Power spectrum ot the fluctuating P«o» e^«' «o "* 5 
HflW MJ/kfl. E -fraction ol total fluctuation energy IQMz). 

nressure signals. and . dominant frequency was observed in ead) 
«7the cower spectra. Figure 5 shows a typical power spectrum of 
1 floating component of Ore pitot signal The power density £ 

is normalised so that 

jcy = i •••<“> 


Tiling 'll* lotiriltunoi bo' 1 ' * idet ol C«l' 1 » l ’ on 42) »nd Uko differ- 

collating yield* 

d^ = dp +2 ^ • (4) 

■- -- 

*j==sr£^-tasscas 

cr-X.'SEi s»sr2— - a-— -• 

dp = d« ...(5) 

where a is the local speed sound By introducing Mach number 
into Ibe above equation yields 

dp _ m— . . (6) 

P = ■ 

Substituting this into Equation (4) gives 

a.(M.r)* • •> , > 

p t V U 

•p* fluctuating components ?\ and * «»" **" j ** “ Sml " 
perturbations ! p , and du- Equation (7) then Kcomes 


Further analysis sbowe’d that this dominant frequency was pro- 
oortional to the square root of the flow temperature (Fig. 6). Susa 
!he soeed of sound in Ibe lest gas was proportional to ihe squire 

Jcngih D. ibe cbaracurisik frequency would be 


a V^RT _ tf ...(U) 

J ’~ D~ D 

where T and R are assumed to be constant. Generally. Equation 
(ID can be written as. 

/,=cVr •••**> 

where the constant of prt>f»rtionality C «corpo^* £ 

° f SLTli?JK?« fld^T^- through 
repression analysis. The origin of the dominant frequency in 
frees tream disturbance has not been investigated but ,t 

Uted that it may be pertinent to the fluctuation of the some surface 

At the call of ibe nozzle ibrosi. n 

tL freestream dUturban.es discussed above were eventually 
-I 


£ l 0000 


< 9 > 

\"V"M42l P, 

, t can thus be seen that for hyP*«^ flow ' W ^ o/pto. 
velocity fixation is U 

ss&E £ 

Fourier transform was performed on me 



Figure 6. Paak frequency as a function of temperature. 
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- FWire 7 Peak mass Dow fluctuation as a tur>c6ono( 

present authors. R « (R # r*- 
received and amplified by -* 

plate model. kadingto ttansiUon “ .tout the 

fiction process, tbe boundary * . defined Reynolds 

frequencies earned in the djs preferred, or tbe most 

number and Mach number. * rt fene rally depends on Macb 
amplified, frequency. 1,1 jugnalion lempera- 

^Vr’^rSV^rdn —» «■ * used to 
lure rabo T. IT* * 7 amplitude A««ltoe neutral 

cakulaie tbe rauo oftbe J ty ^ Uled by tf* turbulent 

point to tbe wlU The catcolatioo is often earned 

boundary Uyer oo toe twule wait toe » ^ vjJuc of tbe 

out in terms of mass flow fluctu , ^ ^ maJ4 flow 

Jocal mass flow fluctuation, u * £ « repfT . 

j fluctuation of tSJ IpUtude ratio A/A* <k ampbfi- 

seots. in a similar manner ' P . vllcou , forcing theory, 

cation of instability waves^ Mac*. * * u fluencies at Macb 
computed tbe peak mas* ** F 7) showed that tbe 

4 5. His numerical result (re pw * ^ bu ^ highest 

most amplified frequency, or ^lo^Reynolds number. As 
amplitude ratio, is relatively bigb^J*^ fluency 

Reynolds number increases. M* m P -rowth. If 

decreases. This trend coincides w,U>.be bound^ ^ 

tbe wavelength of **®?*V^* as tbe boundary layer 

boundary layer iluctneu. frequency or tbe most 

grows with Reynold. velocity of the 

amplified mode deer eases, pro mjder of tbe waves are not 

w.vesrem»imconstmt.N^«b»t ^"^^ „p to the 
distinguished in forcing theory since it »PP<*‘ P 

neutral point. representing tbe am- 

An envelope may be added to tte ct^« ^ ^ p S y ^ 

pbtude ratio for various . . b p^^ular curve repre- 

langenl point of ‘ “^where^that particular mode becomes 
tents tbe Reynolds number where P« . ^ Reyno lds 

tbe most amplified. PJotting r **l fi indicates 

sKaVSsr-cssaiS 


figure 8. The frequency of the most ampkfnd mode as • functon- 
^ Reynolds number. 




It can be seen that the variation of the most amplified frequetcy 

characteristic frequency may be chosen to represent tbe ov» 

frequencies for a large Reynolds number range For 
amplified frequencies tor a g j foury j tu 

boundary layer flow at Macb 3 3, t^oen 

F = 125x10"’. In the present wort tbe numerical result of 

m" ck<‘t was used as a reference even though the nominal Mai 

number of the experiments was sligbUy higher than tbe Mad 

tu ecu*.' ■* *» «*r“i -*r 

r mge covered by the flat plate model. F m was estimated to be > 

‘^mdto^Mtonle as peak frequencies of measured fl^b*** 

. „ i«ife sicnalt, under various test conditions, are plotted ■ 
“Ke dominant frequencies in the freestream noire wen t- 
£*^der of magnimde lower than toe most amphfcd 

“j 

disturbance kve> ,, ». large freestream parameter fluctuation «*•- 
sitiesdo not always mean strong disturbances to tbe bou *0^ 
(r instability and transition. Tbe support of tors condt»M> ^aree 
from Wells's.'* 1 experimen.al investigation of toe 
freestream turbulence™ boundary Uyer transition 
sure gradient at low speed. He compared bs measurements »t> 


. . nf n-nmrtkinality. By differentiating the 

where C is tbe constant of FOP 01 ™""* 1 

.bovx equation toe foUowmg are obtained 
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figure 9. Plot ol the nondmensionalised dominant frequency 
eqainsl stagnation enthalpy. 
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those of Scbuhauer and 

found to* * »** ^^wKpeSed to have differ- 

,1* traiwition Reynolds num , _ 04%, ibe transition 

cot values For f ^ b y 

Reynolds number obl ""^ b * A fcuxot two. To eslim.te 
Scbubsuer and ^^t^Les in «,* d to trsnsi- 
d»e magnitude of the beestr quantity must be ex 

tion, the power spectntt^^ ^ R^ ^ with the 
unined as the wives. ‘ f of iw itioo. In Uct. 

boundary Uy« th,ck ^**, ^ exocrimenully that the major 

S~ -X*» w~.ta.wy ~i« 

tbe flow velocity * t lh ' ^{“^j Jbich travels at the speed or 
frequency tbe wavelenj mucb different from the 

sound relative to to* flow « k fiow . Therefore. a 

wavelength of the TS wave *£ expected, as far as 

correspondent relation Uconcerned. In 

tbe growth of disturbance in ^ of freemeam distur- 

Ibe present investigation, a 1 g P* * ban< j that bad a cone- 
bance energy was disn.butedm a fluency barto 

t pondmg wavelength 0*uc ^ amplification process in 

thickness. It would not be f * model. Disturbances of 

*e boundanr }*** ,° f te deleted maximum frequency might in 
frequencies higher Howt ver the contribution from these 

fact eaist in freestreanv H wouW not he expected 

frequencies .tol to* Act* * dist® bance in free stream was 

srras 

wave number. 


detection of transition 


Many physical quantities 

the boundary layer, and any ^ ^ (]]) the 

criterion for determining ^ s un ton numbers first 

onset of transition as the po However even within 

*» i r" ^£LS Sw h < lm ** *■£ 

not always , fl , p, Jte with uniform pressure, the 

,y. Based on distance from leading edge 
heat transfer rate dec ^“ c V #<)W u diiturbed. the transition 

^mis^fi^Tor^ present wori, as the point where the Sun- 
ton number has a minimum value. 



Fiouw 10 Deviation of measured heat transfer rate from tworefcaf 
v protection at x » 50 mm. 

d,»n tbe expected TS frequencies. This vibration was attributed 10 
!he stress wave travelling through tbe shock tunnel test sectwn 

wall and the model suspension. f ..mentaSv 

Shock lubes by nature create a flow which is fundamenta-y 
unsteady, and care must be taken when applying steady jtuie 
analysis to tbe results of such experiments. Dunng toe penod 
when uncontaminated lest gas is passing over toe (rrodej^ *“£ 
nation and static pressure levels are not constant but decay *i 
time by an amount which is a function of >be mrniel o^mmg 
conditions. However, a region of flow e xuU .bere tbe rroo^f 
pressures across tbe expansion no tile is steady, and toe M*di 
Simiber of the flow can be considered constant. Tbe ratio of sa- 
nation pressure to dynamic pressure will also be constant and f« 
SHTfield around any body shape to berteady. toe locrity 
mcxsmedstatic pressures must be .tody when normalised «to 
respect to stagnation pressure. Provided these pressure ratio* 
eJtoiA, then a small change in the overall pressure level wiD no* 
cause substantial change to tbe overall flow 6dd or the shape of 
toe streamlines, and a steady state analysis msy be userh 
The situation for beat transfer is slightly Afferent. 

U> be considered steady, toe Stanton number must be conrtart 
wbfch implies that tbe beat transfer must decrease in proportion to 
toe product of the mass flow rale with toe recovery eotoalpy. If 
toe tort flow conditions are computed on toe *“ U "'P 0< *' ** 
irentropic expansion follows shock reflection m tbe rtsgnabco 
re*ion. P tben it can be shown that, for constant Stanton number, 
hit transfer must have toe following relationship with stagnation 

pressure P* 


fl “ ^0 *1 


...(16) 


analysis of experimenta l data 

Z*M vta.tta, c ta . C 

taubiuiy -•••■ 1. il Sb.povJ.v~ 

quiet subsonic wind tuorvel * Y correUtio n between 

it was found tost there existed pro v - hrjUonf o( the model, 
the Tollmien-Schhchling wares ^ synchronised with that of 
The vibration frequencies seentod to hV^^chanism which 
instability waves. It wss be ev __ vibrJl -, on , n subsonic bound- 

S^SJSJSS ^ »' taubaity tal 

sition in hypersonic flows. e uansducets mounted 


where y is an effective ratio of specific beats for toe nozzle expm 
Tir^r • 1 of M. beat transfer dependenU on sUgnabon 
pressure to toe power of 1 14 For the purposei of d ‘ lennimn *J b * 
duration of toe steady flow period, this powe, cT A* wasjnored 
for numerical expediency, and toe heat transfer rates were 
normalised by sugnation pressure directly. For a ‘ U *“‘ 

■inn nr es sure of 30% this would lead to an error of 5% in Stamoo 
SE3 ^ substantially affect toe conclusion, drawn 
from the data. Reasonably steady normalised ob L*^ 

in toe laminar flow regions even with this asmplificntjore Fig®* 
10 compares toe measured beat transfer 

from leadine edge with calculated values using Eckert s<» sefa- 
^ tom^ratJl method under different tort condhrons^ 
conditions at toe beginning of the tost period were selectod fa toe 
calculation Tie JLsured beat Uansferratow^ukcna,^ 

moment when toe pressure ratio first reached toe stemt, ^rtate The 
relative errors are seen to be within ±20%. Figure II sbo x a 
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Figure 12. Effect of boundary lay* coding on transition. 
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increase in beat transfer in laminar region due to recombinalionof 
dissoefeted oxygen has been reported for a shap cone model 
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n» (Fit 13) the measured points collapsed closely 00 to a sing 
c^ Vi* UnpUed that, as pointed out by Ro«<*>. to pnmary 

unit Reynold, number coupbng. n both Rg. 12 ar^Fig^ . 
effect of Mach number was not eliminated. «.e. all “* d, “ P° 
were plotted regardless of tbe Macb number, involved in tefl 
flow. The general trend of Fig. 12 show, that transition Reynold, 
number deceased with wall cooling, contra? to to observ.bo. 
made bv Caryt' 1 * but in agreement with tbe result, by 
^triideSTnd Stetson and Rushtoo<'«. Another important 
difference between the present investigation and those of previous 
Si in ibe method of wall «*>liog. In tbi, expenmto 
the variation in wall U> recovery temperanre ratio wa. aebi eved 
by changing the .ugnrtioo enthalpy of tbe flow while the 
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‘T^Ld bv Stetsooand Rusblootw. The method, used by most 
:^K.Ug.«o^«* ^ opposite For this reason. Stetson ard 
Raton's conical da, a were plotted ,n Fig 12 for «mpan^ 
Figure 13 implies that there would be a simple relation t** 
transition Reynolds number and umt Reynolds number, and An 
relation incorporates to effects of Macb number, pressure and 

"IftsUowing NagelV*" discussion, the conelauon between Re, 
»nd Re. can be written as 


Re, = aRe, 


04 


..( 21 ) 


Transition onset Rey^d, numbers. Kmcbon of unit 
Reynolds number. 


where the constant a is determined through tbe experiment asd 
was found to be 2687. The Re, -Re. relationship was compared 
with experimental and theoretical results of other investigator. ■ 
Fig. 14. The transition zone was correlated using Dbawin and 
Nirasimha’rf 16 * proposal, it. 
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ABSTRACT 


In response to the need for ground testing facilities for super 
orbital re -entry research, a small scale facility has been set up a 
the University of Queensland to demonstrate the Superorbital Ex- 
pansion Tube concept This unique device is a free ptstor > dnvtn. 
triple diaphragm, impulse shock facility which uses the enth py 
multiplication mechanism of the unsteady «F*ns.on 
the addition of a secondary shock driver to further heat the dn 
eas. The pilot facility has been operated to produce quasi-steady 
test flows in air with shock velocities in 
It km/s and with a usable test flow duration of the order or 
15 its. An experimental condition produced in the facility with to- 
tal enthalpy of 108 MJ/kg and a total pressure of 335 MPa is 
reported. A simple analytical flow model which accounts for non- 
ideal rupture of the light tertiary diaphragm and the resulting 
entropy increase in the test gas is discussed It is shown that 
equilibrium calculations more accurately model the unsteady ex- 
pansion process than calculations assuming frozen chemistry. This 
STKE the high enthalpy flows pcixfuced in the facihty can 
only be achieved if the chemical energy stored in the test flow 
during shock heating of the test gas is partially returned to the 
flow during the process of unsteady expansion _ Measurements of 
heat transfer rates to a flat plate demonstrate the usability or the 
test flow for aerothermodynamic testing and comparison or these 
rates with empirical calculations confirms the usable accuracy or 
the flow model. 


notation 

a 

speed of sound 

d 

diameter of tube 

h 

specific enthalpy 

l 

slug length 

l m 

asymptotic slug length 

M, 

reflected Mach number 

m c 

m D y 

mass per unit area of gas 

mass per unit area of tertiary diaphragm 

Pr 

Prandil number 

q 

heat transfer rate 

Re 

Rey nolds number 
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temperature 
time 

u velocity 

W shock density ratio 

X, shock lube length 

p Mirels parameter 

y ratio of specific heats 

p. dynamic viscosity 

p density 

Subscripts 


1 . . .20 gas slates 

variation of gas state within flow region 
e edge of boundary layer 

/ interface 

r recovery 

R reflected shock 

$ shock 

H- wall 


1, INTRODUCTION 


With the increasing interest in the design of vehicles intended 
to enter the atmosphere of other planets or return to the Earth’s 
atmosphere from beyond Earth orbit, experimental facilities will 
be required to allow aerothermodynamic ground testing of flight 
vehicle concepts and components. The envelope of flight regimes 
|or these vehicles is expanding. Recent parametric studies on 
missions which involve Mars entry and Earth return 1 " foresee 
aerocapture manoeuvres involving atmospheric approach veloci- 
ties of up to 14 5 km/s. It has not been possible to produce flow s 
suitable for aerodynamic testing, at such speeds, in any existing 
ground based facility. Even the most capable testing facilities such 
as free-piston driven shock tunnels and expansion tubes are limit- 
ed in the maximum flow enthalpies they can produce. In the case 
of reflected shock tunnels this limit is due to the high temperatures 
and pressures which must be contained in the stagnation region 
before expansion of the flow through the nozzle, and the associat- 
ed energy losses through radiative cooling. The standard single 
stage driver expansion tube does not experience this limitation, 
but its performance is limited by the strength of the shock w ave 
produced by its driver. Flow velocities up to 9 9 km/s have been 
achieved in a free-piston driven single-stage driver expansion 
tube' 3 '. 
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Non-reflccted shock tunnels can produce shock speeds which 
are fast enough for superorbital studies*'*. However, all the energy 
is added to the test gas across the shock, which produces a disso- 
ciated. ionised plasma quite unsuitable for aerodynamic testing. 

The Superorbital Expansion Tube is a new hypervelocity facili- 
ty designed to produce flow conditions in a variety of test gases at 
velocities exceeding Earth orbital velocity. It operates as a free 
piston driven, triple diaphragm, impulse windiunnel. It differs 
from the standard expansion lube layout with the addition of a 
secondary driver section in which the primary driver gas is used to 
shock heat a secondary driver gas, boosting the driver perfor- 
mance. A small scale Superorbital Expansion Tube facility is 
currently in operation at the University of Queensland 11 *, and is 
being used to investigate these hypervelocity flow regimes using 
both air as test gas. reported here, and carbon dioxide test gas. 

Test conditions arc only useful if the state of the gas in the test 
flow can be quantified. As well as discussing the operation and 
performance of the facility, this paper also outlines the equilibri- 
um anahsiv used to calculate the ^tate ot the test flows used in 
experiments in the small wale Superorbital Expansion Tube. This 
analysis is also appropriate for flow calculations in a standard ex- 
pansion tube or any other facility using an unsteady expansion 
generated from the shock induced rupture of a light diaphragm 
with a large pressure difference across it. 


2, OPERATION OF THE SUPERORBITAL EXPANSION 
TUBE 

The operation of the Superorbital Expansion Tube is described in- 
Morgan and Stalker*', h uses the phenomenon of enthalpy multi- 
plication of the test flow through an unsteady expansion as used in 
the standard expansion lube, described originally by Resler and 
Blossom 1 *'. In this process the total enthalpy of the flow is deter- 
mined by the extent of the final unsteady expansion of the test gas. 
The layout of the facility is shown in Fig. I. The tunnel is dri- 
ven by a helium filled, free piston driver which is used to shock 
heat a helium filled secondary driver section (2111 m in length), 
separated from it by a heavy diaphragm (1*6 mm mild steel. 
52-5 MPa hurst pressure!. The primary driver, designated TQ. was 
originally designed for use with a small scale shock tunnel. It is of 
limited performance as it can only produce driver gas tempera- 
tures up to I800K. The compression tube of the primary driver 
has an internal diameter of 100 mm and uses a 25 kg free piston. 
The secondary driver, shock tube and acceleration tube sections 
are all of constant area (38 mm internal diameter). 

The facilitv can be operated in two driver modes. In the first, or 
reflected, mode the driver gas is stagnated at the end of the sec- 
ondary driver section and exhausts into the test gas contained in 
the shock tube section (of adjustable length, initially 3 15 m). 
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Figure 2. Wave diagram of the flow process in the Superorbital 
Expansion Tube. 

which is separated from the secondary driver by a thin steel 
secondary diaphragm (0-6 mm, 12-4 MPa burst pressure!- In the 
second or straight-through mode the thin steel secondary di- 
aphragm is replaced with a light cellophane diaphragm (23 mm. 
95 kPa burst pressure) which ruptures with the arrival of the 
shock. The second mode is the preferred mode of operation. Be- 
cause the primary shock speed is overtailored, the reflection of the 
shock in the first mode would be expected to contaminate the 
driver gas and reduce the driver's performance. The experiments 
described in this work use the preferred straight through mode of 
operation. 

In both modes the resulting shock wave propagates into the 
shock lube section and accelerates the quiescent test gas. The 
shock wave traverses the shock tube and ruptures the light tertiary 
diaphragm (9 mm grocery wrap material, 20 kPa burst pressure! 
and then accelerates to a higher velocity as it passes into 
low pressure accelerator gas (helium or air) contained in the 
acceleration lube (1-289 m in length, adjustable). The test gas 
following the shock is also further accelerated as it is expanded 
unsteadily into the acceleration tube. At the exit of the 
acceleration tube the tfesi gas passes over the instrumented mode! 
mounted in the lest section. Helium is the preferred accelerator 
gas as it has a low density and subsequently lower Mach number 
giving lower pressure and temperature ratios than air for the same 
shock speed. 

It is the addition of the secondary driver section which, under 
the appropriate operating conditions, is able to boost the perfor 
manee of the facility beyond that of a standard free piston driven 
expansion tube as reported by Neely er u/ 1 -*. In an expansion tube, 
the performance of the facility, that is the flow enthalpy which can 
be produced, is a function of the strength of the initial incident 
shock wave produced by the driver. Morgan and Stalker 1 *' de- 
scribe the performance enhancing mechanism of the secondary 
driver with the use of a wave diagram of the flow (Fig. 2). 

The primary and secondary driver gases in regions 2 and 3 
respective!} . which are both helium, have the same prewure and 
velocity as they are separated only by a gas interface. Thus if the 
free piston driver is operated so that the primary shock is overtai- 
lored then the helium in region two will have the higher speed of 
sound The driver gas is therefore able to drive a higher shock 
speed through the test gas contained in the adjacent shock tube 
than would be possible for a single driver arrangement where 
region three drives the test gas shock directly 
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Flow Condition 

secondary driver fill pressure 
shock tube fill pressure 
acceleration lube fill pressure 
test gas 
accelerator gas 
shock tube length 
final primary shock speed 
final secondary ***** *f** d 
final tertiary shock speed 



ioeMJ/Kg 

(kPa) 

60 

(Pa) 

800 

(Pa) 

to 


air 


He 

(ml 

1 625 

(m/s) 

6450 

(m/s) 

7280 

(m/s) 

13000 


Table 1 

Experimental condition 


j establishment of flow conditions using air 

TEST GAS 


A number of test conditions for the small scale Superorbital Ex- 
pansion Tube have been established by investigating a range of 
filling pressures in the secondary driver, shock lube and accelera- 
tion lube sections. The highest enthalpy experimental condition 

produced in the facility is set out in Table I . 

For the high shock velocities produced in the shock tube 
(> 8 km/s) it was observed, during initial testing, that there was 
severe attenuation of the flow with the shock being slowed by up 
to 2 km/s when it reached the tertiary diaphragm. The degree of 
attenuation is a function of the shock velocity. To avoid the result- 
ing decrease in performance due to this viscous flow attenuation, 
sections of the shock tube were removed to reduce us length (to 
1625 m) so lhat the secondary shock reaches the Kill ary 
diaphragm before significant attenuation has occurred. Final 
tertiary shock velocities in excess of 13 km/s were achieved, for 
the condition described in Table I. where a penod of usable 
quasi -steady test flow of 15 ps duration was observed (Fig. 5). 

Static pressure levels were recorded at a number of stations 
down the length of the three sections. These measurements were 
also used to calculate the shock speeds in each section. In the ac- 
celeration tube, ionisation deteciors were mounted in parallel with 
the pressure transducers and were used to confirm the shock 
speeds A typical record of wall static pressures measured in the 
acceleration lube is shown in F.g. 3. illustrating the development 
of a region of steady flow between the tertiary shock wave and the 
trailing expansion wave. Pauli '' proposed that it is either this 
expansion wave, or the reflected w ave emanating from the interac- 
tion of the upstream head of the unsteady expansion with the 
driver- test gas interface in the shock tube, which terminates the 

ICM flow. . 

A standard PCB piezoelectric pressure transducer, in the 

configuration described by Neely n at'-', was initially used to 
measure the centreline phot pressure levels in the test flow. 
Shielding was required to protect the expensive PCB transducer 
from fragments of the primary diaphragm which may be earned 
down the tube by the flow after rupture. This shielding increases 



Fioure 3 Acceleration tube wall static pressure histories showing 
y calculated time for the passage of test slug 

(tertiary shock speed - 10-7 km/s). 
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Fioure 4 . Detail of aluminiurrvTead bar gauge used to measure 
centreline pHot pressure. 

the response rime of the gauge, due to the filling of the cavity be- 
tween the transducer and its shield, to the order of 10 ps (com- 
pared with 2 ys for unshielded operation as was observed for the 
flush mounted wall pressure gauges). This configuration is there- 
fore unsuitable for the high enthalpy, xhon duration flows of ini* 
studv. The shielded PCB was used for set up runs to determine il a 
suitable flow was established, in which the bar gauges 
described below might be used. 

The bar gauge has been developed as a relatively cheap dev ice 
which can therefore be placed directly in the flow without shield- 
ing The bar gauges used were based on a design developed at the 
Australian National University* 1 for use in the T3 free-piston 
driven shock tunnel. They consist of a thin piezoelectric film 
sandwiched between a sensing bar of aluminium and a tail of 
tin-lead alloy (Fig. 4). The device operates by measuring the in- 
tensity of the stress wave which propagates down the sensing bar 
to the piezoelectric film when the gauge is aerody namically 
loaded. The material and dimensions of the tail are chosen such 
that the stress wave which propagates into the tail will not reflect 
back dow n its length before the period of test flow is past. 

Typical wall static pressure and centreline pitot pressure histo- 
ries. recorded at the exit of the acceleration tube, are shown in 
Fig. 5 for the high enthalpy lest condition. The sialic pressure 
trace shown is from an upstream transducer as the exit transducer 
trace was found to be more susceptible to noise, as can be seen in 
Fig 3. High frequency oscillations visible on the traces are consis- 
tent with ringing of the gauges in the mounts and were also seen 
in calibration tests. This behaviour is not thought to be a flow 

phenomenon. _ . 

Calculations indicate that the passage of the slug of accelerator 
gas (5 mm in length) should last less than a 5 ps, and is therefore 
not resolved by the instrumentation. The traces arc therefore 
indicative of the passage of the test gas. The calculated pressure 
levels for these test flows are indicated in Fig. 5 and are discussed 



Figure 5. Centreline p»tot pressure history (upper plot) and 
(lower) wall static pressure history for 108 Mikg test conation 
(Note: P t measured 780 mm upstream of P) 
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in * later section As can be seen from the pitot pressure trace, 
problems were experienced with the rise time of the bar gauges 
used to measure the pitot pressure. Calibration tests with the bar 
gauge flush mounted in a blanked off tube indicate a gauge rise 
time of 4 ps. This rise lime is increased with the addition of the 
flow establishment time for a gauge placed in the freestream. 
While this is not ideal, it is adequate for resolving the region of 
steady flow which is observed in the static pressure traces. The bar 
gauges would also be expected to suffer less from temperature 
effects in the hot flow due to the displacement of the piezoelectric 
film from the sensing surface. Development of the bar gauges, 
which are custom made, is continuing in house, and an improved 
response is anticipated, which should give better definition of the 
limits of the test flow. It is also possible that the use of scored 
primary diaphragms will eliminate the problems with debris and 
allow the use of non shielded PCBs to measure the pilot pressure 
levels in the flow . 


4. ANALYSIS OF TEST GAS FLOW STATES 


4.1 Approach 

This section covers the analysis of the processes to which the test 
gas alone is subjected, and leads to an estimate of the flow condi* 
tions which are achieved in the expanded core flow. No attempt is 
made here to model the driver flow although the perfect gas analy- 
sis from Morgan and Stalker 5 ' was used to help in determining 
the optimum driver conditions. The analysis is one -dimensional, 
and not fully predictive in nature, in that the experimentally 
measured values of the secondary and tertiary shock speeds are 
used. This is seen as being an essential tool for the preliminary 
evaluation of the experimental results. 

In a parallel study 4 * 1 an axisymmciric Godunov scheme with 
laminar viscosity and equilibrium chemistry is being used to sim- 
ulate the experiments, and is expected to add important detail to 
defining the extent of the usable flow in space and time. 

The approach of oilier researchers to calculating expansion tube 
conditions is interesting, and can potentially be usefully applied to 
the current configuration to improve understanding of the com- 
plete flow conditions. Wilson 110 ' has performed a one -dimensional 
finite rate chemistry simulation with diaphragm inertia represent- 
ed by a finite holding time for which the shock lube acts as a 
reflected shock tunnel. The code was used to simulate the condi- 
tions produced in Hypulse, the GASL expansion tube facility* 11 ', 
and predicted an unexpectedly high level of dissociation in the test 
gas due to the strength of the reflected shock from the secondary 
diaphragm. Subsequent analysis by Bakos and Morgan* 1 -' with an 
improved diaphragm anal) sis. and a criterion for the onset of 
chemical freezing, indicated that this level of dissociation was 
overestimated. As discussed below, chemical freezing of enthalpy 
has extra significance in the Superorbital Expansion Tube because 
the enthalpy multiplication effect across an unsteady expansion 
does not act on frozen chemical enthalpy. 

Jacobs 4 * performed a perfect gas axisvmmctric analysis of the 
H)pulse facility, with viscous effects, which revealed useful infor- 
mation about the integrity of the contact surface, the extent and 
uniformity of the test flow disturbances, and departures from one 
dimensional flow in general. The high static enthalpies involved in 
the superorbital facility preclude the use of perfect gas analysis, 
though the usefulness of a two-dimensional code is clearly 
illustrated in this report. 

It is seen that none of the currently available codes is able to 
fully predict the operation and test conditions of the superorbital 
flow. The anal) sis developed here was configured to model the 


dominant effects influencing the flow (i.e. viscous shock anemia- 
lion, diaphragm inertia and the unsteady expansion with near 
equilibrium chemistry) separately and sequentially in a way which 
realistically follows a sample of gas through its passage to the test 
section. Noting the limitations of this approach, good agreement 
with the experimental results was observed. 

4.2 Implementation of the flow model 

As stated above, the observed shock speeds in the shock and ac- 
celeration lubes are used in conjunction with the know n quiescent 
test gas pressure to calculate the secondary and tertiary shcci 
strengths. Real gas calculations of these conditions are made using 
an equilibrium analysis. The static pressure and the entropy level 
behind the shock are determined and these are used to determine 
the full gas state via an equilibrium scheme. The validity of 
assuming an equilibrium state is discussed later. 

Calculation of the changing state of the flow through the 
unsteady expansion is implemented in a stepwise manner. 
The expansion process is divided into a large number of velocity 
increments (for these calculations. 100) which proceed from the 
calculated starling velocity of the test gas (state 8, Fig. 2) to the fi- 
nal velocity after expansion (state 12). This final test gas 
velocity is set equal to the observed tertiary shock speed as it can 
be shown by calculation that full Mircls 1141 development of the ac- 
celerator gas slug length occurs due to its low quiescent density 
Calculation across each increment is governed by the standard 
equation for a simple unsteady expansion wave which relates the 
change in static enthalpy to the change in velocity and the local 
speed of sound. 

d/r = -ad# 

A calculation of the entropy and static enthalpy of the flow is 
made at the initial state eight. For equilibrium or frozen flow the 
expansion is Isentropic. The gas condition is uniquely defined by 
the static enthalpy and static entropy and the Equate code 4 *' 1 is 
used to iterate, at the end of each computed step, a solution for 
pressure, temperature and gas composition. The Crek subroutines 
developed by Pratt and Wormcck' 1 ** are used to calculate the equi- 
librium composition at each point. These subroutines are valid for 
35 ionized and unionized species in air between 300K and 24 000 K- 
The process proceeds until the final velocity is reached, thus sup- 
plying the final stale of the test gas. In this process it turns out that 
the determination of the initial gas slate after diaphragm rupture, 
which then undergoes unsteady expansion, is a critical step. 

The lightest possible tertiary diaphragm is used to separate the 
air test gas and the helium accelerator gas. Ideally this diaphragm 
would rupture instantaneously and. having no mass, would have 
no further influence on the flow. But as has been reported in previ- 
ous expansion tube work’ 17 ' even the lightest diaphragm docs not 
behave ideally resulting in two principal non ideal effects. These 
are reflection of the incident secondary shock wave into the on- 
coming test gas and momentum loss from the expanded flow due 
to the energy required to accelerate the diaphragm mass after 
rupture. This behaviour can be observed on a more detailed wave 
diagram for the flow in the vicinity of the tertiary diaphragm (Fig. 6). 

Operation of the facility has show n that the tertiary diaphragm 
material docs not petal, but is sheared off around its circumfer- 
ence. leaving the now unattached material, whether it be solid or 
vapour, to travel with the flow. The subsequent motion and 
behaviour of the diaphragm is unclear. The strength of the reflect- 
ed shock as it traverses the test gas is described well by a simple 
inertial model, but the values of tertiary shock speed measured 
could not be achieved if the diaphragm residue remains as an 
obstacle to the flow for more than about one tube diameter down- 
stream. Stalker 1 * 1 has suggested that fragmentation may allow test 
gas to leak past the diaphragm. Independent of whatever happen** 
to the diaphragm, the boundary conditions for the test gas across 
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Figure 6. Wave diagram for the flow near the tertiary diaphragm. 

ihe unsteady expansion are determined from the strength of the 
reflex-led shock and from the experimentally measured tertiary 
shock speed. The analysis presented here relates these boundary 
conditions to giv e a prediction of the final flow properties. 

Petaling of the diaphragm may be the preferred rupture mecha- 
nism as it would anchor the material to the diaphragm station, but 
for such a light material it is hard to achieve In this instance the 
diaphragm is impulsively loaded to two orders of magnitude 
above its static rupture pressure, and circumferential shearing 
occurs. Attempts to actively rupture the diaphragm in a controlled 
manner have been made- 1,1 and continuing work on these methods 

is being carried out in this department. 

These non-ideal diaphragm effects must be accounted for and 
have a significant influence on the calculation of the test gas state. 
The slug of shock heated lest gas which arrives at the ternary di- 
aphragm (Fig. 6. region 6) is processed by the reflected shock 
wave (region 8) and there is a resulting entropy rise in the test gas 
which depends on the strength of the reverse shock It then passes 
through the unsteady expansion which develops behind the 

accelerator gas-test gas interface. 

The source of shock heated test gas for the expansion is non 
uniform. The first pan of the test gas slug to arrive (region 6a) is 
stagnated by the shock which is reflected from the tertiary di- 
aphragm (region 8a) and is then expanded from this state. Once 
the diaphragm ruptures though, this stagnation of the flow no 
longer occurs as the reflected shock is attenuated by the expansion 
waves (res ions 8 eventually to a Mach wave, and is no 
longer stro'ng enough to bring the oncoming ICM gas to a hall. The 
test gas may thus be expanded from a range of initial gas stales 
z) all of varying entropy. 

Calculation of the expansion process and the resulting state of 
ihe lest gas must be made for this ranee of starting conditions. It is 
perhaps easiest to consider the two hounding extremes of the start- 
ing flow. The first extreme is the initial pan of the test gas slug to 
enter the expansion (region 8u) which will also become the first 
part of the final test flow This region of test gas is taken as initial- 
ly stationary with conditions determined behind a fully reflected 
shock bringing the lot gas to resi The other extreme is the termi- 
nating case (region S/il which is the final pan of the test gas slug 
io make it into region 12 before the end of ihe acceleration lube. 
The calculation of this state is less straightforw ard. 

Bv tracing a particle path through the unsteady expansion, the 
exit point of the particle may be determined. Pauli' - ’ set out some 
simple one step analytical relations for the panicle path through 
an unsteady expansion, directly relating the entrance and exit 




1 Figure 7. Tracing the terminating particle path and the reflection 
from the interaction of the secondary interlace and the reflected 
shock, through the expansion for the 108 MJAg condition. 

trajectories. Unfortunately these relations are for an idea! gas and 
assume a constant value of the ratio of specific heats. As the 
variation of y has significant effect through the unsteady 
expansion, it is necessary to fully track the particle using the equi- 
librium scheme described earlier in this section. By adjusting the 
entrance time into the expansion, different particle paths may be 
simulated, and the appropriate one chosen. 

The disturbance (u + a wave) emanating from the point where 
the interface, separating the test gas and secondary driver gas, 
catches up with the attenuated reflected shock must also be 
tracked as it, rather than the trailing expansion wave, may temu 
naie the test flow (Fig. 7). This occurs in a similar fashion to the 
process described by Paull« 7 > whereby the test flow is terminated 
by the wave emanating from the interaction between the 
driver-test gas interface and the leading expansion wave. 

Calculation of the terminating particle paths, for the flow condi 
tions and facility geometries described here, reveal that only a 
very small fraction of the shock heated test gas in region 8 reaches 
the final lest flow and that it is the trailing expansion wave which 
terminates the test flow rather than the reflected wave (Ftg. 7). 
This indicates an expected period of steady flow of 17 ps which is 
of the order of that observed. This also suggests the tunnel geome- 
try described may not be the optimum configuration for the flow 
condition. Disregarding any further attenuation of the ternary 
shock, calculation of the trajectory of the reflected wave suggests 
that the optimum length for the acceleration tube at this condition 
is 5 m. This would significantly increase ihe slug length of test 
flow and thus extend the lest lime to 60 ps. However, operation of 
the facility with longer Acceleration tubes was found to cause 
severe attenuation of the tertiary shock speed. 

To determine the point of interface catchup and the strength of 
the reflected shock at ihe beginning of the terminating partivle 
path, it is necessary io model the behaviour of ihe Ionian 
diaphragm after rupture. Meyer-'*’ developed an analytical one- 
dimensional model for the impact of a shockwave on a movable 
wall such as a diaphragm, and the subsequent wave behaviour af- 
ter impact. The model though is only accurate for weak incident 
shocks. A more appropriate anal>tical model has been developed 
for the Superorbital Expansion Tube and is described here for the 
strong shocks experienced in the flows considered in this study. 

The residual influence of the tertiary diaphragm after rupture is 
represented by the inertia of the diaphragm as it is convected 
downstream'*’ in a similar form to the action of a free pivton. The 
mass of the ruptured diaphragm is accelerated downstream by the 
pressure difference across the diaphragm, the pressure behind ihe 
reflecied shock being far greaier than the quiescent accelerator 
gas pressure (initially as much as 700 000 limes) and 
the subsequent shock heated accelerator gas pressure. As this oc- 
curs expansion waxes propagate upstream la>m the diaphragm, 
eventually attenuating the reflected shock to a Mach wave. The 
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accelerating diaphragm also sends compression waves down- 
stream into the accelerator gas which quickly converge to form 

the icniary shock (Fig 6). . . . 

As the reflected shock and the diaphragm are initially in i close 
proximity (within one tube dtameier) and the sound speed beh.nd 
the reflected shock is very high, it is assumed that the expansion 
waves emanating from the diaphragm immediately catch up with 
the reflected shock and that the attenuated shock pressure is im- 
mediately realised on the face of the diaphragm. Thus we have 
specified a region of uniform gas velocity between the reflected 
shock and the diaphragm remnants. The velocity of the interface. 

formed by the diaphragm remnants, is determined by the 
integration of its acceleration from rupture. 

, V , ,, Oh-Pii). ...<2> 

' ' (»'Ol +n, G) . 

where m.n is the mass of the diaphragm material per unit area and 
m„ is the mass of the gas per unit area w hich is contained between 
the reverse shock. 5*. and the tertiary sh.vk. 5,. 


= P S (.T|-X s ,) + P"(» f , -'i) 
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The diaphragm velocity may then be used with the known in- 
coming test gas velocity to determine the strength of the reverse 
shock at each moment'' 1 . 
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A calculation of the xtual slug length of lest gas at the end of 
the shock tube is made by solving Equations t6f*' and this is 
compared, for the test flow described in Table I. to the calculated 
asymptotic slug length in Table 2. 

It can be seen from Table 2 that for both conditions the slug of 
test gas does not reach full development in the shock tube, and 
thus the velocity of the interface and the adjacent test gas does not 
reach the attenuated shock velocity. Mirels uses Equation (7) to 
calculate the flow non uniformity between the shock wave and the 
contact surface. 

b ... (7) 


As the reverse shock traverses the test gas it induces an entropy 
chance in the gas which may be determined from the Mach num- 
ber M,. calculated from Equation (4). The period it takes the 
shock to traverse the slug of test gas is small (< 15 (is) and it can 
be shown that only a small displacement of the diaphragm will oc- 
cur in this time indicating that it is permissible to determine the 
reverse shock strength from the theory presented. It should also be 
noted that as the initial pressure loadings on the tertiary di- 
aphragm are small (less than '/:sth of the burst pressure) there will 
be no significant pre-deformation of the diaphragm further sug- 
gesting one dimensional behav iour and circumferential failure or 

the light diaphragm material. , . 

As this simple diaphragm analysis assumes ideal gas behaviour 
there will be some inaccuracy associated with its use on the high 
enthalpy conditions considered here. This analysis though, is on y 
used to determine the shock sirengih at the beginning of the termi- 
nating panicle path. The use of an adjusted value ol the ratio of 
specific heats, which is set equal to that necessary to match the 
initial value of M, calculated from Equation (4). against that re- 
quired by the iterative equilibrium scheme to fully stagnate the 
te-t g.t' i region Rat. helps to minimise the error. It should be noted 
that! he eniropy change across the reflected shock (M* - 3 6) is 
much less than that across the incident shock (M„ - .0). so that 
errors in this procedure arc second order effects. 

To calculate the velocity of the gas behind the attenuated re- 
verse shock at the beginning of the terminating particle path (Stale 
git) i he velocity immediately upstream of the shock must be 
determined Although the shock heated lest gas in region 6 is 
nominally indicated as being a region of uniform flow in rig. - 
and Fic 6. it ■' in fact a region in which there is a velocity gradi- 
ent induced by viscous effects. Mirels' analysis'”’ was used to 
determine the velocity of various points in the region and this is 
found to have a significant effect on the final values or pitot 
pressure calculated for the expanded test flow . 

Mirels' theory implies that due to the viscous effects of the 
boundary laser grow ing behind the incident shock, the vehxity of 
the gas behind the shock will approach the attenuated velocity of 
the shock at an asymptotic slug length determined by Equation <5) 


_P 

(P. 


Bi-.i-W 5 

{>-) 


If we assume that the flow density remains essentially constant, 
an estimate of the velocity at a distance behind the shock, which is 
less than the asymptotic slug length, may be made. 






..( 8 ) 


This is used lo determine the incoming lest gas velocity adja- 
cent to the beginning of the terminating panicle path (region 8n). 
Jacobs ,n ’ has shown numerically, using an axisymmeinc Navier- 
Stokes code, that the rise in velocity behind the shod is essentially 
of the form described by Equation (8). 

For the initial velocity determined in this manner, the gas 
velocity and gas state across the reflected shod is calculated 
(Table 3> and is then expanded through the unsteady expansion 
Once the starting conditions for the expansion process are de- 
termined the remnants of the diaphragm are disregarded and it is 
assumed they have no further effect on the flow sure. The effects 
of the curv ature and shifting of the origin of the expansion fan arc 
essentially confined to the immediate vicinity of the tertiary di- 
aphragm and for ease of calculation are also disregarded. The 
unsteady expansion is assumed to behave as a simple linear fan 
for the implementation of the equilibrium scheme described 
earlier in this section. - 
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asymptotic slug length 
actual dug length 
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Table 2 


1 - — 

Flow condition 

108 MJ/ 

v* 

Flow slate * 

a 

n 

reflected Mach number 

pressure (bP*) 

temperature 

static enthalpy (MJ/Vg) 

velocity <rTV ^ 1 

ratio of vpwific heats 

gas constant (J/VgK) 

total pre^'-ure" tMPa) 

total enthalps • (MJ/Vg) 

366 
7672 
16 620 
582 
0 

1 38 
571 
7672 
582 

3 51 
7095 
15 431 
556 
332 
1 39 
570 
7140 
55 7 

• determine J b> iNcmropic stagnation 

Table 3 

Calculated flow states behind the reflected shock 
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Flow condition 

108MJAg 

calculation method 

equilibrium 

frozen 

flow state 

12 

a 

■ 

a 

velocity T 

(m/s) 

13000 

13000 

,13 000 

pressure 

(Pa) 

6028 


30 

temperature 

<K> 

5976 

Uilifl 

174 

density 

xlO-^kg/m') 

0222 

0218 

0030 

static enthalpy 

(MJ/kgi 

236 

22 8 

289 

Mach number 


657 

6-66 

29-0 

gas constant 

(J/VgKl 

454 

448 

571 

ratio of spectre heats 


143 

M3 

1 64 

pitot pressure 

(LPji 

375 

368 

51 

total pressure* 

(MPa) 

2868 

3354 

— 

dissociation fractions 

N 

065 

062 

1-00 


O 

100 

1-00 

UX) 

'equal to measured tertiary shock speed 




•determined b> i sent topic stagnation 





Table 4 




Final test flow states calculated using equilibrium 


and frozen theory 
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late the limits of the range of test gas flows produced at the end of 
I he unstead) expansion. The calculated gas states of the test flow 
(region 12) at the end of the acceleration tube, are set out in Table 4. 

The resulting calculation of static pressure and pitot pressure 
levels for each starting condition can also be seen in Figs 5a and 
5b. superimposed on those pressures measured in the experimen- 
tal flow. Examination of the flow states before (Table 3) and after 
expansion (Table 4) indicate a decrease in the static pressure of 
the test gas compressed by the reflected shock due to the attenua- 
tion of the shock after diaphragm rupture. A slight downward 
pressure gradient in lime remains after expansion to the final flow 
state. It may be seen from Fig 5 that although the overall levels of 
predicted pressures are close to observed values, the experimental 
static and pitot pressure levels appear to increase slightly through 
the lest flow. 

Although the analysis included viscous effects behind the sec- 
ondary shock, the expansion to the tertiary shock speed was per- 
formed using inviscid theory. It is thought that boundary layer 
growth in the expanded test gas (region 12) may be causing the 
noted discrepancy. A calculation of the displacement thickness of 
a laminar air boundary layer at the end of the test gas slug 
<0 97 mm for the 108 MJ/Vg condition) indicates a contraction of 
the flow core sufficient to produce a net pressure rise in the test 
flow*. The analysis of Akman and Morgan 1 * 1 is targeting this 
specific problem. 

It is interesting to compare the final flow states calculated by 
modelling the expansion as an equilibrium process with calcula- 
tions made assuming a frozen expansion process. Calculations of 
the expanded state of the test gas. required to match the observed 
final velocity, were made for a frozen expansion by fixing the 
chemical state of the lest gas after processing by the reflected 
shock (slate 8). The resulting expanded frozen stales are also set 
out in Table 4. 

The high dissociation fractions of the nitrogen component of 
the air for the frozen cases indicate that large amounts of chemical 
energy are stored in the flow after shock healing in the shock tube 
section. It can be seen that the resulting temperatures, static 
pressures and pitot pressures of the expanded frozen test gas are 
unrealistically low gixing support to the assumption of an approx- 
imately equilibrium process occurring through the unsteady 
expansion To achieve the high enthalpy flow conditions produced 
bv the Superorbital Expansion Tube in experiments, the flow ener- 
gy used to dissociate the test gas when it is shixk heated in the 
shock lube must be at least partially restored to the flow through 
recombination of the test gas through the unsteady expansion. II 


this energy remained frozen as chemical energy the observed test 
conditions could not be achieved. It can be seen that allowing 
the expansion to proceed in an equilibrium manner results in 
significant recombination of the nitrogen component. 

In the Supcrorbitat Expansion Tube, as in all facilities which 
use the unsteady expansion process, the total enthalpy is not fixed 
at any point in the flow as it is in shock tunnel facilities. Rather, 
the Superorbital Expansion Tube relies on the enthalpy multiplica- 
tion mechanism of the unsteady expansion and thus any increase 
in energy due to recombination through the unsteady expansion 
will be multiplied by the expansion to produce the elevated firuJ 
enthalpies recorded. This is one of the main advantages u hich ex- 
pansion tubes have over more conventional shock tube facilities in 
producing high enthalpy flows. 

It should be noted that in reality, the chemical reactions occur- 
ring through the unsteady expansion do so in such short distances 
and limes that finite rate chemical kinetics would come in to play 
This is especially true for the beginning of the lest slug where ihf 
test gas traverses the expansion wave near its origin and thus ex- 
periences the highest rate of expansion. The equilibrium model 
described here, though, allows a straightforward calculation giv- 
ing a useful indication of the final flow states, without having to 
consider the complexities of finite rate chemical kinetics. It i* 
planned that future development of this model will include some 
accounting for these chemical kinetics via the approximate 
scheme described by Bakos and Morgan 1131 . 


5. HEAT TRANSFER RATE MEASUREMENTS TO A 
FLAT PLATE 


As a demonstration of the usability of the superorbital velocity 
flow produced in the facility for aerothcrmodvnamic testing, ex- 
periments in the test flow were conducted to measure heat transfer 
rates to a flat plate for the 108 MJ/kg condition. The flat plate. 
65 mm in length and 25 mm wide, was mounted along the centre- 
line of the acceleration lube exit at zero incidence to the flow (Fig. 
8). It was instrumented with six platinum-quartz thin film resis- 
tance thermometers, flush mounted along the surface at 8 mm 
intervals beginning 13 2 mm from the leading edge. Although the 
duration of the test flow is brief there are sufficient model lengths 
of test gas in the test slug to establish steady flow over the model 
and a typical unsmoothed integrated heat transfer rate history 
illustrating the essentially steady nature of the test flow over the 



Figure 8. Flat plate geometry lor heat transfer rate experiments. 
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Figure 9. Typical inlegrated heat transfer rate history for 108 MJ/Vg 
condition, measured 13-2 mm from the leading edge. 

As local Reynolds numbers for the flow over the flat plate were 
below 5 x 10 6 . calculations of the expected heal transfer rates for 
a laminar boundary layer were made using Eckert's empirical 
reference enthalpy method'? 1 *. 

-- .- 1 (9) 

*. =0 332p«,Pr 'Re 2 (h,-h m ) 

The starred values are evaluated at a reference enthalpy defined 
by. 

*• = 0 S{h, + h.)+0 22(ft, + »,) 

The calculated test flow state, set out in Tabic 4, was used to 
calculate conditions at the edge of the boundary layer and at the 
surface of the flat plate using a pressure distribution determined 
by the weak interaction principle for a cold wall 1 --*, which is ap- 
propriate for the geometry and flow state in question. These in 
turn were used to calculate the heat transfer rate to the flat plate. 
The calculated rates are compared in Fig 10 with the surface tem- 
perature histories which were integrated, in the manner of Schultz 
and Jones'?' 1 , to obtain the heat transfer rates. 

For laminar boundary layers with similar profiles, the coeffi- 
cients in Equation (9) are all constant except for the local 
Reynolds number. In Fig. 10 the Stanton number is normalised by 
which is proportional to ^/Rc* 4 **, and the steady level when 
plotted against the wetted length is seen to be constant, indicating 
laminar flow. The level is also seen to be in agreement with 
Ecken's empirical correlation. This is an interesting result as the 
correlations are unproven at these flow velocities and would not 
be expected to account for high enthalpy heat transfer mecha- 
nisms such as radiative healing which increase in significance at 
the elevated temperatures encountered here 

Dips in the surface temperature traces at the arrival of 
the shock, as can be seen in Fig. 9. were evident on many of the 



Figure 10 Plot of St normalised by r 4 * * versus distance from lead- 
ing edge for 108 MJ/Vg 


experimental traces. He? 4 ' observed similar dips in his flat plate 
experiments and postulates that they may be due to a decrease in 
resistance across the gauge due to conduction through the ionised 
boundary layer flow. For a uniform lest flow the initial resistance 
change due to this conduction is large with the starting of the sur- 
face flow and results in a decrease in the recorded voltage. The 
surface flow quickly becomes steady and the subsequent voltage 
increase is due only to the increasing surface temperature. He sug- 
gests that a simple first order correction may be made by effec- 
tively shifting the origin of the heat transfer rate trace lo the base 
of the dip. This would increase the measured healing rates to a 
level greater than the laminar calculation predicts. 

The correlation of the empirical calculations and the integrated 
measurements is good, indicating that the method of determining 
the final test flow states described above is effective and that a 
useful core of lest gas at superorbital velocity has been created in 
the facility. 


6. CONCLUSIONS 

The combination of an added secondary driver section to further 
shock heal the driver gas and the enthalpy multiplication effect of 
the unsteady expansion of the test gas allow the unique Super- 
orbital Expansion Tube facility to produce flows at superorbital 
velocities, without the dissociation problems experienced by con- 
ventional shock tubes at these speeds. The small scale facility 
used to investigate the concept has been operated to produce 
quasi-steady test flows with shock velocities in excess of 13 km/s 
and with a duration of usable test flow of approximately 15 
Such high enthalpy flows are of importance in the design of any 
vehicle used to operate within a planetary atmosphere at velocities 
exceeding Earth orbital velocity, such as in proposed missions 
reluming to Earth from Mars. No other existing ground based 
test facility can produce test flows at these velocities, with gas 
composition suitable for aerodynamic testing. 

Consideration of the effects of boundary layer grow th and non 
ideal diaphragm rupture on the flow have allowed a better under- 
standing of the dominant flow mechanisms and enabled more 
accurate matching of calculated and observed test flow states. In 
particular account has been taken of the increase in the entropy of 
the test gas associated with the shock reflection at the tertiary di- 
aphragm and the resulting non uniformity of the test gas slug in 
the shock tube. 

It has also been shown that much better agreement u ith experi- 
ment is achieved if the unsteady expansion is modelled as an 
equilibrium process rather than as a frozen chemistry process, 
as the high enthalpies produced by the facility are not possible un- 
less some of the chemical energy stored in the flow ufter shock 
heating of the test gas is relumed to the flow during the unsteady 
expansion process. 

Heat transfer rate measurements to a flat plate at zero incidence 
in the superorbital flow were of the same level as those predicted 
by Eckert's empirical correlations for laminar flow. These results 
show that the lest flow produced by the facility can be used for 
meaningful aerothermodynamic research and further confirm that 
the analytical method of determining the final test flow states is 
both reasonable and useful. 

These initial results have demonstraied the viability of the 
Superorbital Expansion Tube concept and while further testing on 
the small scale facility will continue, work is now proceeding on 
the design of larger scale facilities It is foreseen that such facili- 
ties will be invaluable for the hasic aerothermodynamic research 
required for future interplanetary missions involving atmospheric 
manoeuvring as well as for the design and testing of vehicles to be 
used on such missions. j 
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Balances for the Measurement of Multiple Components of Force in 

Flows of a Millisecond Duration j eje^ j /pcfpQ 

... a ™ i«n c 1 Tuttle 3 and J.M. Simmons 4 
It M^J^l EngUrme. The University of AUSTRALIA 

Abstract Tt.il tapS“ il l n °" s 

components of force - lift. <“S p ^ d ^ basics G f the design of the 
have a duration of about conc sc t at incidence in 

balance are presented and resultso .. clions> These results indicate that the 

the T4 shock tunnel are com P^ . 9 L 220 mm long model. Also presented are 
prototype balance performs appSon of .he rfeconvolutiou force 

s u rrhe" h n. 

Key words: Force measurement. shock lunnel. iustn.menta.ion 

1. Introduction 
ne measurement of forces o„ 

resiricled by the short durations fo^hj^ desiening balances for use in 

representative flow. Ptogres ! h»* ^ ^ (J( J e „ Gronig 1993; Naumann et 

nows of duration as short as a “ for measuring a single force component 

at- 1993; c arbonruo 993 y Anewwq ^ recently a. the Umversity of 

the deconvoluuon drag balana. n 11iis technique involves intetpreung 

Queensland (Sanderson and S * 5ting conmc led to the test model to 

transient signals from sbamgg -phe deconvolution drag balance has 

“ - Sh0rt - 1 - 

(Porter et al. 1993). , . r orcc balance is extended to the simultaneous 

,n “ ffl# The basie design of this balance is 
measurement of hit, arag an v * balance for measurements with a 15 

considered and the instal at ‘ on ® f (hJs modc i at sma ll angles of incidence in a 
semi-angle cone is discussed. R twL A further application of the 

. »— * — * 

also discussed. 

2. The three-component balance design 

* rW revolution force balance consists of a single, 2 m long sting 
The three-component deconvolutio Each of the short bars is 

" — 2&2X* a. iu L point. A suain gauge bridge is 
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also attached to the sting 200 mm from the modelling junction. Combinations of the 
u ain sTanals ftom the four bars are used to produce two output signals - one 
reuionding primarily to a lift input signal and the other responding primarily to a 
pitching moment input signal. The strain measurement in the sting responds i primarily 
m a drag input. Inevitably there is some coupling amongst these output signals. 



pig. 1 Schematic of the three-component balance arrangement. 

The time histories of the three outputs related to lift, drag and moment. y L {t), y D (t) 
and v^r) can be related to the time-histones of the lift force, drag force and pitching 
moment on the model. «.</>. u„«> and via nine impulse response funct.ons. Ttas 
coupled convolution problem can be written in matrix notation as in Mee et al. ( I ) 


f % 

yt 

/ 

y D 

s 

O’*#- 

V 


Gld Glm 
G dl G dd G Dt4 

’ML GmD G m m 


« l 
u t 


I 


where the y vectors are formed from the discretized output signals and the u vectors 
^e formed form the applied load time histories. The square G matrices are formed 
from the impulse response functions. G # being the impulse response for the y t output to 
a u, input. If there is no coupling amongst the outputs then the off-diagonal 
submatrices in the impulse response matrix will be null. Further details on the des.gn 

of the force balance are given in Mee et al. (1992). . 

The nine impulse response matrices were found from a senes of bench tests in 
whi ™a weight was attached to the model by a fine wire. The wire was cut. thus 
producing a ftep change in the load applied to the model. The output signals resulung 
from such a load change were recorded and processed to produce the impulse 
responses The linearity of the system enables the responses to flow-type loa mg 
distributions to be determined by superposition of the results of several tests for single 

loads applied at various locations on the model. .... 

In experiments in the shock tunnel, each of they, outputs was measured and tune 

domain, coupled deconvolution techniques were used to determine the time histones 
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„f the lift and drag forces and filching moment on the model. The experimentally 
££& *£ response functions were used for this deconvolutton. 

3. Tests of a prototype balance 

. . was installed in a 220 mm long, 15° semi-angle, aluminium 

Tte prototype balan«w« tn^alled^ ^ ^ ^ ^ ^ devicc was fi rs , tested 

cone as indicated in Fig. • hlng wcight s the model. Results for a single 

^"d‘,0 L conewere^aredjith £ ^^",ng 

Unbalance was "stlnsiuve to the distribution of load on the model. 



Fig. 2 Details of the connection of the balance 
to the sting. 


Fig. 3 Typical Pilot pressure trace for present 
tests. 


rwsrfnrrrwvl in the T4 free-piston shock tunnel (Stalker and 

EXpe i988Ho test toe^rformance of the balance. The sting was supported by fine 
Morgan 1988 ^ ^ in thc planc D f the lift force, and was 

How so that ihe only acrodynanuc ^ 

~ ..... t( ,, ts t u e cone was set at three angles of attack - 0.0 , 2.3 ana a.u . 

“"'fm the tests nitrogen, at Ihe conditions given in Table 1. was used A typtcrj 
, Pit0 , 0 „ ssure during a shot of the tunnel is given in Fig. 3. Note tba 
filtered trace ot rttot p & , . Qr v,: tr9n/ The Pitot trace shows a small 

- mean level 

""'ft 'SZlmX Pi™ wavesTniltd by the arrival of thetes. gas 
It takes aboul tsuu p P reflect from , he free end and return to the strain 
lo traverse the g , .. . fl ow has traversed approximately 15 

STe^tt C°de convolution of .he mcasuted sigoals continues beyond ibis point 
model 1 g ■ .1 resoonses at these larger times deteriorates. 

but the accuracies Jp „r»o“ a geTnals obtained from Ute ftve strain gauge bridges 
Examples of the r^ voliage^g^ ^ ^ ^ ^ ^ show „ for ^ 

are shown in tg- • transverse bars joining the model to 

S,"r«S measured in the sfutg itself. The signaHcnoise ratio is 
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fru.nrt to be adeauale using semiconductor strain gauges. The relatively large 
oscillations in til signals in the transverse bars are associated with tbe dynamtes of the 
model/sling configuration. 


i hiAfnmsi conditions * test cas isnitrogen. — — — . 
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Deconvolved signals are shown in Fig. 5 for the three incidences of the present 
tests The deconvofution process leads to amplification of noise on signals and the 
results ^hown have been passed through an 8-pole Butterworth filter with a cut-off 
r no of fi kHz The drag signals show an initial overshoot during the nozzle 
Sn time hoover t£ fifths, which have a slower response, do not show this 
overshoot Both lift and drag show reasonably steady values during the test time of 
500 800 ns after flow start. The moment signal for the model is defined to be zero 
f . ,f) t forcc acts through a point one-third of the cone axial length from the base 
^"e cone mZSSL option of the lift force for an inviscid conical fl<£ 
The results indicate that the centre of pressure during the test time is about 5 mm (2% 

.he computed data oflo« 

(1969) There is a small difference in the ratio of specific heats (1.4 in Jones 1969 and 
( 1 32?n the experiments) but a comparison of the results, shown in Table 2 gives an 
ndk fion of the ^performance of the balance. There are some uncertainties in the 

conditions of the flow in t>* « 

draglre'about' 10% lower in the experiments but the ratios of lift to drag agree well. 
At This stage the accuracy of the balance has not been quantified. 


5 Thrust measurement on a 2D scramjet nozzle . . 

m deconvolution force balance it also being used to measure thmst produced by a 
symmetrically divergent nozzle with 1 1“ ramp walls. The area ratio is approximately 


no 






i riinoc 5n which the measurement of tensile 
- S wS - «W thrust. The system is show, in Fig. 6. 



rig. ^ - 


Test No. 


3745 

3752 

3748 


Incidence 

(deg.) 


0.0 

2.5 

5.0 


Experiment 
Drag. D Lift, L 170 
(N) (N1 


510 

460 

520 


20 

210 

420 


0.04 

0.46 

0.81 


Computation 

Drag. D Lift, L L/D 

(N) gQ — 

5 » 0 0M 


530 

590 


220 

470 


0.42 

0.80 


. . . ono nM , t0 t hc nozzle to provide rigidity and 

The stings are twisted through 90 shown |his to ^ acceptable. There 

bending stiffness. Finite element ") . ch st j ng . The static pressure along the 

are strain gauges before and after thc plane will be measured in the 

samp wans and the £-£££ by the pressures and dta. 

experiments. Any diff halance will be due to skin friction, 

determined with the dec ° nv0,ut ' 0 " )a J “ |at on for a tunnel-type loading have been 
The strains from a finite c,cr ^"^ response based on another finite element 
deconvolved successful y . m ^ u factured and instrumented and benchtesting 

computation. The model has hce comoutations (Fig. 7). Differences between 
shows encouraging agreement wi P attributed to different end 

the step responses after arrival of Based on ^se results and 

^“^«eS«o„volu.ion balances i. is ettpected that .he present balance 

will perform satisfactorily in tests in T4. 

6. Conclusions 

o iKr^-comoonent deconvolution force balance have 
The results of tests using a satisfactorily in a shock tunnel flow of a 

indicated that the balance p erf domci we || on a model that was not 

millisecond ti^del was 220 mm long and had a mass of 1.9 kg. 

S£fS5f liw >° «-* - ° r,he ba, ““- 
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Fig. 6 Schematic of the thrust balance. 
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Fig. 7 Experimental and computed step 
responses for a point load on the nozzle. 
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